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SLOW CONVERGENCE IN PERIODIC HOMOGENIZATION 
PROBLEMS FOR DIVERGENCE TYPE ELLIPTIC OPERATORS 

HAYK ALEKSANYAN 


Abstract. We introduce a new constructive method for establishing lower bounds on 
convergence rates of periodic homogenization problems associated with divergence type 
elliptic operators. The construction is applied in two settings. First, we show that solu¬ 
tions to boundary layer problems for divergence type elliptic equations set in halfspaces 
and with infinitely smooth data, may converge to their corresponding boundary layer 
tails as slow as one wish depending on the position of the hyperplane. Second, we con¬ 
struct a Dirichlet problem for divergence type elliptic operators set in a bounded domain, 
and with all data being C°“-smooth, for which the boundary value homogenization holds 
with arbitrarily slow speed. 


1. Introduction 


The focus of the paper, as the title suggests, is on quantitative theory of periodic 
homogenization of divergence type elliptic operators. Lately, there has been much interest 
around effective estimates on convergence rates for homogenization problems associated 
with linear elliptic operators in divergence form, see for example [11]-[12], [2]-[4], [17], [16], 
[6]. A particular common feature of these papers, is that they all establish upper bounds 
on the speed of convergence of homogenization, in other words these papers, among other 
results, measure how fast the homogenization holds. However, results showing limitations 
of the speed of the process, i.e. estimating to which extent homogenization may decelerate, 
given that the homogenization takes place of course, seem to be extremely scarce in the 
literature. A few instances of this type of set-up around the divergence setting which 
we are aware of are the following. It is shown in [3]- [4], that the Dirichlet boundary 
value homogenization in U’ can not be faster than a certain algebraic rate depending on 
1 < p < oo and on the geometry of the domain (see [3, Theorem 1.6] and [4, Theorem 
1.3]). The next one studied in [19] is related to boundary layer problems set in halfspaces, 
and shows that depending on the position of the halfspace, convergence of the solution 
to its boundary layer tail can be slower than any algebraic rate (see [19, Theorem 1.3]). 
Finally, in [8] it is proved an existence of one-dimensional examples in almost periodic 
homogenization with fixed boundary and source terms and oscillating coefficients, where 
homogenization of solutions in L? is not faster than a polynomial rate. 

Here we will be interested in developing tools that will address how slow the convergence 
can actually be. We start with the discussion of the hrst problem considered in this paper, 
and will introduce part of the key ideas in that setting. 

For a scalar a G M and a unit vector n G consider the following Dirichlet problem 


-V • A{y)Vv{y) = 0, y G VLn,a := {y G : y • n > o}, 
v{y) = VQ{y), y G dVLn,a- 


The main assumptions concerning (1.1) which will be in force throughout are: 
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• (Periodicity) The coefficients A and boundary data vq are Z'^-periodic, that is 
for any y and any /i G Z*^ 

A{y + h) = A{y) and vo{y + h) = vo{y), 

• (Ellipticity) There exists a constant A > 0 such that for any x G and any 
y G one has 

A|xp < A{y)x < A“^|xp. 

• (Regularity) All elements of A and boundary data vq are infinitely smooth. 

We will refer to (1.1) as boundary layer problem. These type of problems emerge in 
the theory of periodic homogenization of Dirichlet problem for divergence type elliptic 
operators with periodically oscillating coefficients and boundary data. Understanding the 
well-posedness of problems of the form (1.1) in a suitable class of solutions, and more 
importantly the asymptotics of solutions far away from the boundaries are one of the 
key steps toward obtaining quantitative results for homogenization of the mentioned class 
of Dirichlet problems. For detailed discussions concerning the role of (1.1) we refer the 
reader to [11], [12] and [19]. Below we will briefly review some known results concerning 
boundary layer problems. Interestingly, the asymptotic analysis of (1.1) depends on certain 
number-theoretic properties of the normal vector n. 

Rational directions. We say that n G is a rational vector, and write n G 
if n is a scalar multiple of a vector with all components being rational numbers. One 
can easily see that if n G has length one, then n is rational iff n = ^/|^| for some 
non-zero ^ G Z'^. In this case it is well-known (see e.g. [5, Lemma 4.4]) that there exists 
a smooth variational solution v to (1.1), which is unique given some decay conditions on 
the gradient, and such that there is a constant for which v{y) —?• exponentially 
fast as y • n —>• oo, where the convergence is uniform with respect to tangential directions. 
Although having these nice convergence properties, the drawback of rational directions is 
that the constant may actually depend on a, i.e. translating the hyperplane in the 
direction of n may lead to different limits at infinity. 

Diophantine directions. Following [12] for a unit vector n G set to be the 
operator of orthogonal projection on the hyperplane orthogonal to n. Fix Z > 0 so that 
(d — 1)Z > 1 and for k > 0 define 

(1.2) = {n G : |P„x(OI > for all ( € \ {0}}. 

A vector n G is called Diophantine if it is from for some k > 0. Clearly elements 
of Ak are non-rational directions. Also, it is not hard to verify (see [12, Section 2]) that 
(t(S'^“^ \ Ak) < where o is the surface measure of the unit sphere, and C is a 

constant depending on 1. Thus, the last inequality shows that almost all directions are 
Diophantine. 

Behaviour of (1.1) in the case when n is Diophantine has been studied only recently 
in [11], where it was proved (Propositions 4 and 5 of [11]) that there exists a smooth 
variational solution v to (1.1) which is unique, given some growth conditions, and such 
that for some constant v°° one has v{y) —>■ v°° as y • n —)• oo. Here convergence is locally 
uniform with respect to tangential directions, and is faster than any polynomial rate in 
y • n. Moreover, the effective constant v°° (the boundary layer tail) depends on direction 
n only, and is independent of a in contrast to the rational case. 

Non-rational directions in general. Here we consider directions which are irra¬ 
tional, i.e. from the complement of Observe that not all irrational directions are 

Diophantine, therefore the previous two cases do not cover the set of all possible 
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directions. In a recent work [19], it was proved that (1.1) has a unique smooth varia¬ 
tional solution satisfying certain growth conditions, for which one has convergence toward 
its boundary layer tail far away from the boundaries (see Theorem 2.1 for the precise 
statement). However, the result of [19] does not provide any estimates on the rate of 
convergence given this generality on the normals. It does however show that for irrational 
directions which are non-Diophantine (meaning they fail to satisfy (1.2) for any choice of 
parameters k and I involved in the dehnition), one may have convergence slower than any 
power rate in y • n. More precisely, for a smooth and Z^-periodic function uq : —>■ M 
consider the following boundary value problem 

(1.3) Av = 0 in and u = uq on (9H„, 

where Hn = {x G ; x • n > 0}. Clearly, this problem is of type (1.1) with the matrix of 
coefficients equal to 2 x 2-identity matrix. Then [19, Theorem 1.3] shows that if n ^ 
is arbitrary non-Diophantine direction, then for any p > 0 and any R > 0 there exists a 
smooth function uq : —>■ M and a sequence Xk oo such that if v solves (1.3) with 

boundary data vq, then for each A: = 1, 2,... and all y' G d^ln H H(0, R) one has 

(1.4) \v{y + Xkn) - v°°\ > X'^^, 

where the constant v°° is the corresponding boundary layer tail. Let us note that the left- 
hand side of (1.4) converges to zero as /c —>■ oo, since as we have just said, for irrational 
directions solutions converge to their boundary layer tails. The proof of (1.4) constructs vo 
with Fourier spectrum supported in a subset of I? on which the normal n fails to satisfy the 
Diophantine condition. Then choosing coefficients having an appropriate decay, combined 
with the special structure of the spectrum of vq, immediately leads to the conclusion. We 
stress that this slow convergence result of [19] works for any irrational non-Diophantine 
direction, however, it leaves out the question whether one can go beyond algebraic rates 
of convergence, and perhaps more intriguing, the case of Laplace operator does not give 
an insight into the case of variable coefficient operators, since in the Laplacian setting one 
has an explicit form of solutions which essentially determines the analysis. 

Throughout the paper we use the following notation and conventions. 

T'^ unit torus of M'^, i.e. the quotient space M^/Z'^, where Z*^ is the integer lattice, 

unit sphere of 

0{d) group oi d X d orthogonal matrices, 

M'^ transpose of a matrix M 
X ■ y inner product of x, y G 
|x| Euclidean length of x G 

^n,a halfspace {x G : x ■ n > a}, where n G and a G M, 
halfspace Xlnfl, 

Br{x) or B[x, r) both stand for an open ball with center at x G and radius r > 0, 
(s compact inclusion for sets. 

Positive generic constants are denoted by various letters C, Ci, c,..., and if not specihed, 
may vary from formula to formula. For two quantities x, y we write x < y for the inequality 
X < Cy with an absolute constant (7, and x x y for double inequality Cix < y < C 2 X 
with absolute constants Ci,C 2 - 

Throughout the text the word “smooth” unless otherwise specified, means differentiable 
of class C°°. The term “modulus of continuity” is everywhere understood in accordance 
with Definition 1.1. The phrase “boundary layer tail” refers to the constant determined 
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by Theorem 2.1. Domain is an open and connected subset of Euclidean space. We also 
adopt the summation convention of repeated indices (Einstein summation convention). 

1.1. Main results. The first class of problems we will study in this article, is motivated 
by the results discussed above, and the importance of boundary layer problems in periodic 
homogenization of Dirichlet problem. Most notably, we will show that in the case of irra¬ 
tional non-Diophantine normals convergence of solutions to (1.1) towards their boundary 
layer tails can be arbitrarily slow. Next, in the second part of the paper, we will apply our 
methods developed for the analysis of (1.1) combined with some ideas from our papers 
[2]-[4] written in collaboration with H. Shahgholian, and P. Sjolin, to construct a Dirichlet 
problem for elliptic operators in divergence form set in a bounded domain, where boundary 
value homogenization holds with a speed slower than any given rate in advance. 

We now proceed to formulations of the main results. In order to measure a speed of 
convergence we consider the following class of functions. 

Definition 1.1. We say that a function oj is a modulus of continuity if it has the 
following properties: 

• w : [0, oo) —)• (0, oo) is continuous, 

• OJ is one-to-one, decreasing, and lim oj{t) = 0. 

t^OO 

We will at places abuse the notation and instead of [0,oo) may take [co,oo) for some 
Co > 0 as a domain of definition for modulus of continuity. 

For our first result we will impose a structural restriction on coefficients A. Namely 
we assume that 

(1.5) there exists 1 < 7 < d such that dy^A'^°' = 0. 

In other words we require one of the columns of A to be divergence free as a vector field. 
This assumption is technical and is due to our proof of Theorem 1.1. It is used to control 
the contribution of boundary layer correctors in the asymptotics of boundary layer tails 
(see, in particular, inequality (4.34)). 

The following is our main result concerning the slow convergence phenomenon of 
boundary layers. 

Theorem 1.1. Let 00 he any modulus of continuity and let R > 0 be fixed. Then, there 
exists a unit vector n ^ a smooth function uq : —)> M, and a sequence of positive 

numbers growing to infinity, such that if v solves (1.1) under condition (1.5) 

on the operator, and with n and vq as specified here, then for any k = 1,2,... and any 
y' G dfln,o n -6(0, R) one has 

(1.6) |u(y'-kAfcn)-u“| > a;(Afc), 
where v°° is the corresponding boundary layer tail. 

Remark 1.2. Observe that v°° being the boundary layer tail, implies that the left-hand 
side of (1.6) decays as k ^ 00 , and hence the lower bound is non-trivial. Next, notice that 
we have fixed the halfspace on the direction n G by setting a = 0 in Theorem 1.1. 

This does not lessen the generality, since the case of arbitrary a G M can be recovered by 
a change of variables. However, we do not know if in general one can take the sequence 
{Afc} and boundary data vq independently of a. 

Finally, let us note that the result of Theorem 1.1 shows that there is no lower bound 
for the speed of convergence on the set of irrational directions, in other words convergence 
can be in fact arbitrarily slow. This case is in strong contrast with the case of Diophantine 
normals, where one has convergence faster than any power rate. 
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Our next concern is the question of Dirichlet boundary value homogenization for di¬ 
vergence type elliptic operators in bounded domains. Assume we are given a coefficient 
matrix A = : X —)> defined on some domain X C R‘^ (d > 2) and 

having these properties: 

(Al) for each 1 < a, /3 < d we have G C°°(X), 

(A2) there exist constants 0 < A < A < oo such that 

A|e|" < < A|eP, Vx G X, VC G M'". 

For a function g G C'°°(T'^), and a bounded subdomain D <Z X with boundary consider 
the following problem 

(1.7) — V ■ A{x)Vue{x) = 0 in D and Ue{x) = g{x/e) on dD, 

where e > 0 is a small parameter. Along with (1.7) consider the corresponding homoge¬ 
nized problem which reads 

(1.8) — V • A(x)Vuo(a:) = 0 in iA and uq{x) = / g{y)dy on dD. 

Jrd 

Let us emphasize that we do not impose any structural restrictions on A, nor we 
assume that A is necessarily periodic. This is in view of the fact that there is no interior 
homogenization taking place in (1.7). 

Theorem 1.3. Let A be as above defined on X and satisfy (Al)-(A2), and let eo he any 
modulus of continuity. Then, there exist bounded, non-empty convex domains D <Z X and 
D' D with boundaries, and a real-valued function g G C'“(T'^) such that if Ue is 
the solution to (1.7) for e > 0, and uq to that of (1.8), then for some sequence of positive 
numbers strictly decreasing to 0, one has the following: 

a) \ue^{x) - uo{x)\ > uj{l/ek), Vx G H', A: = l,2,..., 

b) |ue(x) — uo(x)| —5-0 as e —)■ 0, Vx G iA. 

This result should be compared with [2]- [4], where under certain geometric conditions 
on boundary of the domain (such as strict convexity, or flat pieces having Diophantine 
normals as considered in [4] ) it is proved that periodic homogenization of boundary value 
problems for elliptic operators in divergence form holds pointwise, as well as in norm 
where 1 < p < oo, with an algebraic rate in £. However here, we see that again due to 
the geometry of the domain, convergence can slow down arbitrarily. Thus relying merely 
on the smoothness of the data involved in the problem, one can not get a meaningful 
quantitative theory for homogenization problems with divergence structure as above. 

2. Preliminaries on solutions to boundary layer problems 

The aim of this section is to give a precise meaning to a solution of problem (1.1). 
The well-posedness of (1.1) in non-rational setting was first studied by Gerard-Varet and 
Masmoudi [11], and later by the same authors in [12], and by Prange in [19], all in 
connection with homogenization of Dirichlet problem. Here for the exposition we will 
follow mainly [12] and [19]. Let us also note that the results presented in this section hold 
for strictly elliptic systems, however we will only use them for scalar equations, and thus 
formulate the results in the setting of a single equation only. 

Keeping the notation of problem (1.1), fix some matrix M G 0(d) such that Me^ = n. 
Then in (1.1) make a change of variables hy y = Mz and transform the problem to 


( 2 . 1 ) 


-Vz ■ B{Mz)Vzv{z) = 0, Zd>a, 
w{z) = vo{Mz), Zd = a, 
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where v(z) = v{Mz) and the new matrix B is defined hy B = M'^AM. From the definition 
of M we have M = where Ai is a matrix of size d x (d — 1). Then the solution to 

(2.1) is being sought of the form 

v{z) = V{Nz',Zd), where V{-,t) is Z'^-periodic for any t > a, 
and z = {z', Zd) G x M. This leads to the following problem 


( 2 . 2 ) 




N^Ve 

dt 


■ B{9 + tn) 


y{e,t) = vo{9 + an), 


N^Ve 

dt 


V(9, t) = 0, t > a, 
t = a. 


The authors of [12] then show that (2.2) has a smooth solution V in the infinite cylinder 
T'^ X [a, oo) satisfying certain energy estimates. Moreover, if V solves (2.2), then one can 
easily see that v(z) = V{Nz', Zd) gives a solution to (2.1). The proof of this well-posedness 
result is not hard, but what is rather involved is the analysis of asymptotics of V{-,t) as 
t —)> oo. A proper understanding of this problem for V gives the behaviour of solutions to 
(1.1) far away from the boundary of the corresponding halfspace. In this regard, it was 
proved in [11] that if n is Diophantine in a sense of (1.2), then there exists a constant v°° 
depending on n and independent of a, such that \v{y) — < Caiy • for any a > 0 

as y • n —>■ oo, and convergence is locally uniform with respect to tangential variables. 

Shortly after [11] and [12], a refined analysis of well-posedness of problems of type (1.1) 
was given by Prange [19]. In particular he established the following result. 

Theorem 2.1. (Prange [19, Theorem 1.2]) Assume n ^ Then 

1. there exists a unique solution v G C°°{Qn,a) H L°°{Qn,a) of (1.1) such that 

l|Vn||Loo({y.„>i})0, ast^oo, 


/ <oo, 

J a 

2. and a boundary layer tail v°° G M independent of a such that 

v{y) —)> as y ■ n ^ oo, where y G ^n,a, 

and convergence is locally uniform with respect to tangential directions. 

Let us fix here that by boundary layer tail we always mean the constant to which solu¬ 
tions of problem (1.1) converge away from the boundary of the corresponding halfspace. 


3. Arbitrarily slow convergence for Laplacian 

The objective of the present section is to prove Theorem 1.1 for Laplacian. The reason 
for separating the case of Laplace operator is twofold. First, we will introduce part of the 
key ideas that will be used in the general case of variable coefficient operators. Second, the 
setting of Laplacian is essentially self-contained, and is more transparent in comparison 
to the general case which is based on a different approach. We prove the following result. 

Theorem 3.1. Theorem 1.1 holds when the operator in (1.1) is the Laplacian. 

The proof of Theorem 3.1 is based on a series of observations and preliminary state¬ 
ments. We will use the connection of problem (1.3) with the corresponding problem (2.2) 
set on a cylinder x [0, oo) (cf. [19, Theorem 7.1]). For that fix a matrix M G 0{d) such 
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that Mcd = n, clearly M is of the form M = [iV|n] where is a matrix of size dx{d— 1). 
Then writing (2.2) when the original operator in (1.1) is Laplacian and o = 0 we get 

2 

l/(0,t) = O, t>0, 

tmO)=Po( 0 ), 

where as before V{-,t) is Z'^-periodic for all f > 0, and the action of the operator on V is 
understood as {N^'V 0 ,dt) • {N'^V 0 ,dt)V. As we have discussed in Section 2 , the unique 
solution V of (1.3) (in a sense of Theorem 2.1) is given by 

(3.2) v{y) = v{Mz) = V{Nz',Zd), where y = Mz with y G 0,n and z G 

In this setting the solution V of (3.1) can be computed explicitly, namely we have 

(3.3) V{9,t) = ^ 

where c^(uo) is the ^-th Fourier coefficient of vq. In view of (3.3) it is clear the boundary 
layer tail equals co(uo). We will first establish a slow convergence result for V using which 
we will prove Theorem 3.1. 

Observe that by (3.3) and Parseval’s identity, for any t > 0 we have 

(3.4) \\V{ 6 ,t) — co(uo)|> \\V{9,t) — co(uo)||^ 2 (']pd) = 

^ := §{t-vo). 

CezAfo} 

Proposition 3.2. For any modulus of continuity uj there exists a unit vector n ^ 
a smooth function uq : > M, and a sequence of positive numbers tk oo, k = 1 , 2 ,... 

such that 

§{tk;vo) >u){tk), A: = 1,2,..., 

where § is defined by (S.f). 

As we can observe from (3.4) convergence properties of S depend on the quantity | 
which is the subject of the next result. 

Lemma 3.3. Given any modulus of continuity lv, there exists a unit vector n ^ and 

an infinite set A C Z'^ \ {0} such that for any matrix M = [A'|n] G 0(d) one has 

llV^el <a;(|e|), V^G A. 


(3.1) 


N^V0 

dt 


As usual here as well N is a d x (d — l)-matrix formed from the hrst (d — 1) columns 
of M. Obviously Meu = n. This lemma is one of the key statements used in the general 
case as well. 


Proof of Lemma 3.3. Set := ei = (1, 0,..., 0) G and let Ti C ^ be an open 
neighbourhood of centred at and with diameter less than 

Due to the density of rational directions^ there exists a non-zero G Z'^ such that 
> 2 and 


0 < 


^( 2 ) 




< diam(ri) < 


10 |{( 1)|2 ' 


^For any non-zero v € Q'* the intersection of the ray starting at 0 G R'* and passing throngh v, with the 
sphere is a rational vector of unit length. Hence the density of rational directions. 
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Using the same idea, we then inductively construct a sequence C Z'^\ {0} satisfying 

k < and such that for unit vectors we get 

(3-5) 0< |rfc+i-rfc| < 

for each fc = 1, 2,... . It is clear by (3.5) that for k large enough one has \rk+i — rk\ < 10“^, 
therefore the sequence {rj.} is Cauchy, hence it is convergent. By n we denote its limit, 
which is obviously a unit vector. We claim that n ^ which is due to fast convergence 
of the sequence^ Indeed, assume for contradiction that n G MQ'^. As \n\ = 1, using 

the rationality assumption it is easy to see that there exists a non-zero ^ such that 
n = Co/|?o|- By monotonicity of and ui, along with (3.5) we get 


|(n • eif - (rfc • eif\ < 2|n • ei - r*, • ei| < 2|n - rk\ < 2^ Ir^+i 

j=k 


-rj\ < 






9 10^ 


k = l,2,... . 


By rationality of n we have n • ei = ||^ with po € Z, and for we have ■ ei 
with some pfc G Z. Hence, from the last inequality we get 


Pk 

|^(fe)| 


- pI\U 




Since the left-hand side is an integer and is less than 1 by absolute value for k large enough, 
it must be 0. From here we conclude that the sequence \rk • ei| is eventually constant. By 
our notation this implies 


(3.6) 


Pk ^ , Pk+i 
|^(fc)| “ |^(fc+i)|’ 


k > ko, 


where pk is an integer, representing the hrst coordinate of and ko is a large integer. 
Now, if we have equality in the last expression with minus sign, we get 

1,1, , 2\pk\ 

]_Qfc|^(fc)|2 — — kfc+l ' — l^k ■ ei\ — ) 

which implies that pk = 0, and hence Pk+i is 0 as well by (3.6). We thus see that (3.6) in 
either case of the signs, forces equality within the first components of r^. and r^+i- But 
in this case the same argument with ei replaced by the remaining vectors of the standard 
basis of would lead to equality for all corresponding components of and r^+i. The 
latter contradicts the fact that ^ rk+i which we have from (3.5). Hence the proof that 
n is not a rational direction is complete. 

We now set A = : A: = 1, 2,...} and proceed to the proof of the claimed estimate 

of the lemma. By orthogonality of M for any G A we have 

|^|2 = |M^^|2 = |iV^eP + |n-^P, 


which, combined with Cauchy-Schwarz, implies 

(3.7) = ICI' - In • eP = (lei + In • ^1) (|^| - |n • ^1) < 2|^|2 (l - . 


This is in analogy with a standard fact in Diophantine approximation theory, that the sum of a very fast 
converging series of rationals is irrational. 
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Now choose A: G N such that ^ . We get 

^{k) 


^ . n = • n = • 

Hence by (3.5) we have 


|^(fc)| 






|{<‘)| + ■ (r,+i - rj). 

j=k 


|J(») |£W|| < £|{l*>||ry+i -r,| < 




j=k 


j=k 


10i|^(i)|2 - 2 |^W| 


We thus get 


1 - 


n 


^(A:)| 


|^(fc)| 


< 


Cu2(|e(^)|) 


2|^W|2 

From here, getting back to (3.7) we obtain 
completing the proof of the lemma. 


□ 


The following remarks will be used later on. 

Remark 3.4. One may easily observe from the proof and the density of rational directions, 
that given any r > 1 it is possible to construct A such that for any ^,r/ G A i/ |^| < \r]\ 
then t|^| < |? 7 |. 

Remark 3.5. The set of normal directions satisfying Lemma 3.3 is dense on how¬ 
ever it necessarily has measure 0 if u decreases faster than any polynomial rate. Density 
follows from the proof of the lemma, as there one may start the construction in the neigh¬ 
bourhood of any rational direction on instead of ei. The measure zero claim is due 
to the fact that the set of Diophantine directions, in a sense of (1.2), has full measure on 
and if u! decreases sufficiently fast, then any Diophantine direction fails to satisfy 
Lemma 3.3. 

We now give a proof of Proposition 3.2 based on the previous lemma. 


Proof of Proposition 3.2. Define a;i(t) : [l,oo) —IR+ as follows 


Wl(t) 


1 


t > 1 . 


Here uj~^ stands for the inverse function of w, which exists since u is one-to-one. Moreover, 
without loss of generality we will assume that oji is well-dehned for t > 1, since otherwise 
we will just replace the lower bound of t by a sufficiently large number. It is easy to see 
from the definition of a;i(t) that it is continuous, decreases to 0 as t —)• oo and is one-to- 
one. We now apply Lemma 3.3 for this choice of wi as a modulus of continuity, and let 
n be the normal, and A be the index set given by Lemma 3.3. We dehne vq G C'°°(T'^) 
as follows. First arrange elements of A in increasing order of their norms, i.e. we let 
A = : k = 1 , 2 ,...}, where for k = 1 , 2 ,... . Observe that > k 

for all A G N due to the construction of Lemma 3.3. For ^ G we let be the ^-th 
Fourier coefficient of vq. Next, if ^ G A is the k-th. element of A according to the mentioned 
arrangement, set c±^(uo) = |?|~^, otherwise, if ^ A let c^(r;o) = 0. The sequence 
of coefficients constructed in this way defines a smooth function vq, since the Fourier 
coefficients of vq decay faster than any polynomial rate. Furthermore, as c_^ = G M, vq 
is a real-valued function. 
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By (3.4) we get 


Ht;vo) = ^ 


\2k 


For each A; G N, choose from the condition that 




(3.8) 


A; = 1,2, 

e 


By construction we have tk oo. To prove that §{tk',vo) > uj{tk) it is enough to show 
that > uj{tk), while for this one it is sufficient to prove 

(3.9) > oj{tk), 

since by Lemma 3.3. We now use (3.8) and the dehnition of wi, by 

which (3.9) is equivalent to 


1 > 471011(1^^1)4 = 47ro;i(|C('=)|)o;-' 


471- 


47101-1 ( l|^(fc)|-2|««| 


^-1 Q |^( fc )|- 2 fc ^ . 


But the last expression is equivalent to 


oi 


-1 Q|^W|-2|?«|^ >^-1 Q|^(A:)|-2fc^ ^ 


which holds true, since ^ oi i is decreasing. 

The proof of the proposition is complete. 


□ 


Remark 3.6. It is clear from the proof of Proposition 3.2 that given any 6 > 0 in advance, 
we may drop some finite number of initial terms from A C Z'^, the Fourier spectrum of 
vq, ensuring that ^ for all C £ 

Remark 3.7. By the same way as in Remark 3.5 we may argue that the set of normals 
with the property as discussed in Proposition 3.2, is dense on the unit sphere, however 
with measure zero if uj has a sufficiently slow decay at infinity. It is also clear that any 
prescribed lower bound on the sequence for non-zero f by a given function of 

1^1 would transform to a certain upper bound on the speed of convergence ofV to its tail. 

Proof of Theorem 3.1. Let a unit vector n ^ a Z'^-periodic function vq and 

a sequence of positive numbers {Afc}^i be obtained by applying Proposition 3.2 for the 
modulus of continuity oi given in Theorem 3.1. Also let V be defined by (3.3) for this 
choice of vq. As we have seen in (3.2) the unique solution v to problem (1.3) is given by 
riu) = V{Nz', Zd), where as usual z G Mz = y with M = [A^|n] G 0{d). By 

orthogonality of M we have 

y ■ n = Mz ■ n = z ■ M'^n = z ■ Cd = Zd. 

Thus if we let y = y' + {y-n)n, with y' G cAlru then Nz' = y' for the tangential component. 
We now need to derive some bounds on the new tangential variable z'. Observe that N has 
rank d — 1, hence d — 1 of its rows are linearly independent. Set N' to be a (d — 1) x (d — 1) 
matrix formed by these d — 1 rows of N. From the overdetermined linear system Hz' = y' 
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we have N'z' = y", where y" G ^ is the corresponding part of y'. Using the assumption 
that |y'| < R, we get the following bound 

(3.10) |z'| = |(iv')“^/l < CN\y\ < cnR. 

Now if A C is the Fourier spectrum of vq, by Remark 3.6 we may assume that < 

l/( 8 cAri?), for any ^ G A, from which and (3.10) one gets 

The last expression shows that cos(27r^ • Nz') > V2j2 for all ^ G A and any z' satisfying 

(3.10) . By construction of Proposition 3.2, A is symmetric with respect to the origin, also 
c^(i;o) = C-^(uo) for any ^ G A, and all non-zero coefficients of V are positive and do not 
exceed 1. Hence for any y' G dQn H B(0, R) we get 

(3.11) v{y' + {y ■ n)n) = V{Nz', za) = = 

C6A 

y2 cos{27t^ ■ Nz') > -^§{zd-,vo), 

where S is defined by (3.4). Finally, recall that y ■ n = z^, and hence in the last inequal¬ 
ity restricting to the sequence {tk}^=i, and using the estimate of Proposition 3.2 we 
complete the proof of Theorem 3.1. □ 

It follows from the proof of Theorem 3.1 that we may have local uniformity for slow 
convergence with respect to tangential directions, meaning that the sequence on which 
the convergence is slow, once chosen for the modulus of continuity uj, can be used for any 
R > 0. The only difference is that in this case one should start at a very large index 
(depending on R) in the sequence. 

4. Variable coefficients 

In this section we prove Theorem 1.1 for coefficient matrix A satisfying (1.5). Recall, 
that the previous section established the slow convergence phenomenon for Laplace op¬ 
erator. The main (and big) point that makes Laplacian special in the analysis is that 
one may write the solution to reduced problems (3.1) explicitly. However, for variable 
coefficient case one does not possess explicit forms for the solutions, which necessitates a 
rather different approach. 

We start with some preliminary set-up. For coefficient matrix A by A* denote the 
coefficient matrix of the adjoint operator, i.e. For 1 < 7 < d we let x*’'^ be 

the smooth solution to the following cell-problem 

f-V, • A*{y)VyX*’^{y) = dy^A*^-\ y G 

\/Td x*’'^{y)dy = 0. 

Next, by Vn'^ we denote the solution, in a sense of Theorem 2.1, to the boundary layer 
problem 

. f-V-A(y)Vur(y) = 0, y G 

\vn''^{y) = -x*’^(y), y e dQn- 

Finally we recall the notion of the Green’s kernel. For a coefficient matrix A and a halfspace 
H C the Green’s kernel G = G{y, y) corresponding to the operator —V • A(y)V in 
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domain is a function satisfying the following elliptic equation 


-Vy ■ A{y)VyG{y, ^ = 6{y - y), y G 
<^(^,^ = 0 , yedn, 


for any y € where 6 is the Dirac distribution. To have a quick reference to this situation, 
we will say that G is the Green’s kernel for the pair (j 4,D). Note, that the definition of 
G does not require A to be periodic. The existence and uniqueness of the Green’s kernel 
for divergence type elliptic systems in halfspaces is proved in [15, Theorem 5.4] for d > 3, 
and for 2-dimensional case in [10, Theorem 2.21]. Here we will only use the case of scalar 
equations. It is also shown that if G* is the Green’s kernel for the pair {A*,Q) then one 
has the symmetry relation 


(4.4) G{y,y) = G*{y,y), y,y £ Q with y ^ y. 

From here we see that Green’s kernel has zero boundary values with respect to both 
variables. 

For a unit vector n ^ MQ'’* and a smooth function vq let v be the solution to (1.1) in a 
sense of Theorem 2.1. Set X := y ■ n for y G and let M G 0{d) satisfy Me^ = n. Then 
for any 0 < k < l/{2d) we have 


(4.5) v{y) 


d 2 ,aG'^{n, Mz) X 


A^^{XMz)np+ 


A^'^iXMz) ( 5 j,^x*’“(AMz) + dypvTi^Mz)) 


vo{XMz)da{z) + 0{X-^), 


where G” is the Green’s kernel for the pair and A^ denotes the matrix of coef¬ 

ficients of the homogenized operator corresponding to —Vy • A{y)Vy. Also 52,a denotes 
differentiation with respect to the a-th coordinate of the second variable of G”, and the 
error term 0{X~‘^) is locally uniform in tangential variable y' := y — (y ■ n)n, and is in¬ 
dependent of the matrix M. The asymptotic formula (4.5) is proved by Prange in [19, 
Section 6] for systems of equations. Here, since we are working with scalar equations, we 
have a slightly simpler form of it. 

We are going to switch from the differentiation in y to z-variable. Since y = Mz and 
M is orthogonal, it is easy to see that 

(4.6) Vy = MV,. 

Set ■= X*’°‘{Mz) and ^n°‘{z) := Vn'^{Mz). On the boundary of for each 

1 < a < d we have x*’°^ + = 0 by (4.2), hence taking into account the relation (4.6) 

and the fact that M has n as its last column, we obtain 

+ Vn^){y) = npd,^{x*’^ + 

Since the Green’s kernel has zero boundary data, using (4.6) we have that if 3 y = Mz, 
with 2: G dR'^ then 

(4.7) dy^G^{n,y) = d 2 ,aG^{n,Mz) = n„5,,G°’”(erf, 2), 

where G^’” is the Green’s kernel for the pair MJ^). Observe that G®’” depends 

on the matrix M, however for the sake of notation we suppress this dependence in the 
formulation. For now, dependence of G^’"' on M plays no role, but later on we will need 
to make a specific choice of orthogonal matrices M. 
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Applying these observations in (4.5) and grouping similar terms we obtain 


(4.8) v{y)= d2,dG^’'^{ed,iz',0)) X n^A{XM{z',0))nx 


1 + n„(a^^x*’“(A(z',o)) + a^^v);’"(A(z',o))) 


X vo{XM{z',0))dz' + 0{X-^). 


The next lemma concerns a particular class of integrals of type (4.8). 

Lemma 4.1. (see [1, Lemma 2.3]) Assume H : ^ M. is smooth lA-periodic function, 

n ^ is a unit vector and M G 0{d) satisfies Med = n. Then for h{z') := H{M{z',0)), 
where z! G and any F G L^(]R'^“^) one has 


lim f F{z')h{Xz')dz' = cq{H) f F{x)dx, 

A^oo J^d-l jRd-l 

where cq{H) = fjaH{y)dy. 

This lemma is proved in [1] for functions admitting a certain type of expansion into 
series of exponentials. To obtain the current version, one can take the matrix T in Lemma 
2.3 of [1] to be the identity. 

In the next statement we collect the necessary information concerning the Green’s 
kernel involved in (4.8). 

Lemma 4.2. For any n G and any M G 0(4) satisfying Med = n let G^^{z,T) be 
the Green’s kernel for the pair (M’^A^M, M!^). Set fn,M{z') '■= 52,dG^"'(ed, {z'where 
Then 


(i) fn,M ^ I/r^-i/ n,M(^')^^^1 > 0; 

(ii) sup ||/„,m||l1(iR‘'-i) < oo sup / \fn,Miz')\dfi ^ 0 as A ^ oo, 

n,M n,M|2;/|>y4 

(hi) fn,M £ C^(M'’*“^) and sup ||liURd-i) < oo, where V' is the gradient in 

n,M ’ ^ 


Proof. Recall that G^ is the Green’s kernel for the pair (A^,fl„). The following bound 
is proved in [12, Lemma 2.5] 


(4.9) 


\G^{y,^\<G 


{y • n){y • ra) 

\y - 


y fiy in 


where the constant G is independent of n. It is easy to observe (see e.g. [1, Claim 3.1]) 
that for any M G 0(4) with Mcd = n we have G^f^{z,'z) = G’^{M'^z, M'^F) iov z fiz in 
From here and (4.9), along with the orthogonality of M one has 


(4.10) 


JG' 


0,71/ -^1 

M = 


G^{M^z,M^F)\ < G 


{M'^z ■ n){M'^z ■ n) 
\MTz-M^F\<i 


= G-, 


ZdZd 


z — z 


d’ 


for all z z in with constant G as in (4.9); in particular G is independent of n and 
M. Since G®’"' has zero data on from (4.10) one easily infers that 


\fn,M{z')\ < 


G 

cd - (^;',0)1'^’ 


Vz' G 


which shows that fn,M £ ^) as well as part (ii). For the second statement of (i) 

let P^{z,F) be the Poisson kernel for the pair (M'^A^M,M^). Then for z G and 
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z G we have 


pO,n 


(z,z) 


-(Me,fA°(Me,)d2,,G%^(z, I) 

-n^A°nd^^G°^(z,z). 


The last expression, combined with the fact that^ fg^d {ed,^da{T) = 1, along with 

the ellipticity of AS' completes the proof of the second claim of (i). 

Finally, for (iii) observe that since Gj^ solves an elliptic equation with constant coeffi¬ 
cients, by standard elliptic regularity we have fn,M £ For the growth estimate 

by [20, V.4.2 Satz 3] one has 

|V'{SuOSi“(ej, (/,0)))| < Vi' € 

where G is independent of the unit vector n and the orthogonal matrix M. 

The proof of the lemma is complete. □ 


We now turn to the discussion of the core of averaging process of (4.8). Our next result 
is one of the key lemmas of the current paper. 

Lemma 4.3. Let H be any non-empty set of indices and assume we are given a family of 
functions T = with the following properties: 

(a) F G for any F G T and inf F{x)dx\ > 0, 

(b) sup ||F||pi(j|d-i) < oo and 
F& 


sup / \F{x)\dx —)• 0 as —)• oo. 

F&3^J\x\>A 


Let also u be any modulus of continuity and Sq C he any open subset of the sphere. 
Then, there exists an irrational vector n G Sq, an unbounded, and strictly increasing 
sequence of positive numbers such that for any F G T there is a function vq G 

C°°(T'^) satisfying 


(4.11) 


F{x)vo{XkM{x,0))dx — co{vo) / F{x)dx 


>a;(Afc), A: = 1,2, 


whenever M G 0(d) and Mcd = n. 

Assume in addition to (a) and (b) that T also satisfies 

(c) F G Ci(M'^-i) for any F G T and 

sup ||VF||iimd-l) < OO, 
Fe5 

then the function vq too can he chosen independently of F. 


^It is proved in [12, Sectin 2.2] (see also [19, Section 3.2]) that the variational solution to (1.1) has an 
integral representation by Poisson’s kernel. Moreover, as long as the asymptotics of the solutions to (1.1) 
far away from the boundary of the hyperplane is not concerned, there are no restrictions imposed on the 
normal direction. Since identical 1 is a solution to (1.1) for a boundary data identically equal to 1, we may 
represent this solution by Poisson’s kernel, which shows that the integral of Poisson’s kernel is 1. 






SLOW CONVERGENCE PHENOMENON 


15 


Remark 4.4. Let us remark that the left-hand side of (4-kl) decays as k ^ oo in view 
of Lemma f.l, therefore the lower bound of the current lemma is non-trivial. The Lemma 
shows that under conditions (a) and (b) only, the direction, and the sequence along which 
convergence is slow can be chosen uniformly for the entire family T. Moreover, as will be 
seen from the proof of Lemma f.S, for any F and G from T their corresponding functions 
vo{F) and vo{G) have equal up to the sign Fourier coefficients. 

Proof of Lemma 4.3. For the sake of clarity we divide the proof into few steps. 


Step 1. Construction of n and {Afc}^^. We start by determining a suitable modulus 
of continuity for which we will apply Lemma 3.3 to get the normal n. 

For F G 3" let Ji? be the absolute value of the integral of F over and set tq := 

inf 3p. By (a) and (b) we have 0 < tq < oo. It is easy to see using the fact that tq > 0 

and condition (b) that there exists j4o > 0 large enough such that 


(4.12) 



> 2 / \F{x)\dx + 
J\x\>Ao ^ 


for any F £ S'. For this choice of denote ep := n-Mi- 

F £ T. Since tq > 0 from (4.12) and condition (b) we have 


I\x\<Ao ^(^^dx 


, where 


(4.13) 


0 < en := inf ep < 1. 
FeJ 


We now fix some small constant <5o > 0 such that 


(4.14) I cos(t) — 1| < eo/4 for any t G M with \t\ < Sq. 

Assume that n £ and let M G 0(d) be so that Med = n. We then have M = [N\n] 
where N is the matrix formed from the first {d — l)-columns of M. Observe that for any 
^ G and any x £ one has ^ • M(x,0) = x'^N'^f^. Therefore if for some A > 0 we 

have 27rAAo|A'^^| < do then (4.12) and (4.14) imply 


(4.15) 


F{x) cospvrA^ • M{x, 0)]dx 


> 


' |a:|<Ao 


F{x)dx 


'\x\<Aq 


|F(x)| X I cospvrAx^A''^.^] — l|dx — f \F{x)\dx > -Jp, 

J\x\>Ao 8 


for any F £ T. Define 


uii{t) := 


dp _ 1 

27rAo w-i ilro^) ’ 


t > 1, 


where uj~^ stands for the inverse function of io. Obviously oji is one-to-one, continuous, 
and decreases to 0 as t —>■ oo. It is also clear that ui is well-defined for large enough t, thus 
without loss of generality we will assume that uoi is defined for all t > 1. Applying Lemma 
3.3 for cji as a modulus of continuity we obtain A C and a unit vector n ^ such 
that if M G 0(d) is any matrix satisfying Med = n, then 


<a;i(|C|), V^G A, 


where, as is customary, d x (d — 1) matrix N is formed from the first (d — l)-columns 
of M. Following Remark 3.5 we may assume that n G Sq. We arrange elements of A in 
increasing order of their norms, thus A = : k = 1,2, ...}, where by construction we 

have k < \ for any fc > 1. Moreover, according to Remark 3.4 we may also 

assume that for any A: G N we have 


( 4 . 16 ) 


|^(0|<^|^(A:+1)| 
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where 0 < ^) < 1 is a fixed parameter satisfying 

(4.17) 2sup ||F||ii(]Rd-i)—^ < -^tq. 

Note that the supremum here is finite by (b) and is non-zero by (a). Set 

(4.18) Xk := k = l,2,.... 

It is clear that Xk is unbounded and is strictly increasing. Observe that n, and the sequence 
{Xk} are uniform for the entire family 3". 

Step 2. Construction of vq for fixed F € S'. We proceed to construction of the 
function vq G C'°°(T'^) for the given F G 3“, for which it is enough to construct the 
sequence of Fourier coefficients of vq, which we will denote by {cg(uo)}^gZ‘*- 

Let F G 3“ be fixed. For ^ G if we have ^ G A then set c^(uo) = c_^(uo) = ek{F)\^\~^ , 
where A: G N is the index of ^ in A according to the increasing rearrangement made above, 
and Ek{F) G { — 1,1} will be chosen below. It is important to note that this sign is the same 
for and c_^. Otherwise, if ^ A we let c^(uo) = 0. Clearly the sequence {c^} decays 
faster than any polynomial rate in |^|, hence vq is smooth. Also, since cg(uo) = c_^(i;o) 
for any ^ G Z*^ we have that uq is real-valued. Observe that co(uo) = 0 by construction, 
and expanding vq into Fourier series we get 


(4.19) 



F{x)vo{XM{x, 0))dx 


E 

m=l 


2em{F) 


F{x)cos {2 'kXx'^N '^:= 


2gfc(F) 

|^(A:)|fc 


^k(^) + Zl(A) -|- S2(A), 


where A: > 1, Jfc(A) := f^d-i F{x)cos{27rXx"’"N"’"^^^'>)dx, Ei(A) contains the part of sum 
where m < k and S2(A) respectively sums over the range m > k. In view of the construc¬ 
tion the sums Sj(A), i = 1,2 are real-valued for any A. By the definition of Xk, the fact 
that > k and that u is decreasing we easily see that 27rAA;Ao|A'^^^^^| < do for any k, 
hence applying (4.15) we obtain 


(4.20) 


2 

|^(fc)|A: 


|3fc(Afc)| > 


On the other hand, by (4.16) and (4.17) we easily get 


(4.21) 


|E2(A)| <2||F||^i(r.-i) £ 


m=k-\-l 


|^(m)|r 


< 


3 1^ 

16 


for any A > 1. We now estimate the contribution of the range m < k. The triangle 
inequality implies 


Zi(Afc) -|- 


1^(01 




+ 


Si(Afc)- 






> 


|^(0|- 


:|4fc(Afc)|, 


hence at least one of the terms in left-hand side of the last inequality is not less than half 
of the right-hand side. Taking this into account, we choose the sign Sk in order to get the 
largest term from the left-hand side of the above inequality. This choice of £k, combined 
with estimates (4.20) and (4.21), and the definition of Xk given by (4.18) yields 


2£k{F) 

|^(fc)|fc 


^k{^k) + Zi(Afc) -|- E2(Afc) 


> ^ > 3 1 


w(Afc), 


(4.22) 
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for any k = 1,2,... . The estimate (4.11) of the lemma obviously follows from the last 
inequality and (4.19). 

Step 3. Uniform choice of vq. Lastly, we turn to the proof of possibility of a uniform 
choice of vq under additional condition (c). There is no loss of generality to assume that 

(4.23) —)• oo as t —>■ oo, 

since otherwise we would simply replace a; by a new modulus of continuity ui, where 
oj{t) > oj{t) for all t > 1 and a; satisfies (4.23), by that getting even a slower convergence. 
Thus we will take (4.23) for granted. For fixed F G T, G \ {0}, and A > 0 set 

(4.24) /(A;e):=/ 


Let l<A:<d — Ibe such that the fc-th component of the vector is the largest by 
absolute value. This choice implies |(F'’^^) • ek\ > {d — 1)“^/^|A''^^|, where ek is the fc-th 
vector of the standard basis of Integrating by parts in /(A; in the direction of ek 

we see that 

(4.25) l^(A;OI < 


where the supremum is finite due to the assumption (c). 

By (4.18) we have |^^^^|*'a;(Afc) = |ro for each fc G N. Also, since \k is increasing and 
unbounded by construction, from (4.23) we get XkUj{Xk) —>■ oo as fc —>■ oo. Hence 


(4.26) 


Afc 

|^(A;)|A: 


=: Ufc —>■ OO as /c —)• OO. 


We now choose an increasing sequence of integers {ik)^i where = 1 and if for k > 1, 
ik-i is chosen, we use (4.26) and take ik > ik-i so large in order to get 


(4.27) 


1 ^/d^ 


k-l 




IT 


SUP||VF|L.(.^-)T |^(i„)|.„|A,T;(i.,.)| S 


Clearly the choice of the sequence (ik) is independent of a particular F since constants in 

(4.27) are uniform for the entire family T. We define vq through its Fourier coefficients as 

follows. If for ^ G Z'^ we have ^ for some /c G N then define c^{vq) = C-^(vo) = 

1^1“*'=, otherwise, set cg(uo) = 0. We have that vq is uniform for all F G T. Observe also, 
that Vo is simply the function from Step 2 with the difference that its Fourier spectrum 
is now supported on the frequencies non-zero Fourier coefficients are 

positive. 

We now complete the proof by showing that n, {Aij.}^r vq defined above satisfy 
the Proposition. Plugging vq into (4.19), for each integer A; > 1 we get 

(4.28) / F(x)vo(XM(x,0))dx = j^\(X) + Ei(X) + E 2 (X), 
where Tj., Si and S2 are defined as in (4.19). We have 


k-l 

Si(A) = 

m=l 


1 


|^T7 


/(A; ?(*-)) +/(A;, 
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with I defined from (4.24). From this we obtain 


(4.29) 


(4.25) 

Si(A4)| < 


1 


A 




TT 


k-1 


sup ||VF||^i(]Rd-i) 

F& 


E 

m=l 


1 1 


(4.27) 

< 


1 3 (4-26) 3 1 

< Ye 


For the subsequence 


of A the analogue of estimate (4.21) becomes 


1 oo ^ 

(4.30) |S2(A)| < 2||F||ii(Rd-i) |£(i„)|i™ - 2||^llLi(Rd-i) ^ uim)\m - 

m=k+l m=ik+i ^ 

3 1 < ^ ^ 

16 |^(A+i~i)IA+i~i ^ 16 |^(A)|A 

where as before, we have used (4.16) and (4.17). Finally, the lower bound on (4.28) of the 
Lemma follows by replacing k with in (4.20) and applying estimates (4.29) and (4.30) 
to (4.28). 

The proof is now complete. □ 


Looking ahead let us remark here, that the importance of uniformity of the choices in 
Lemma 4.3 will prove crucial in the applications. We now include a small modification of 
the previous lemma to allow compactly supported functions, as well as shift of the origin 
in the function vq. This situation emerges in applications of Lemma 4.3 to integrals arising 
from Poisson kernels corresponding to bounded domains. 


Lemma 4.5. Keeping the notation of Lemma f.S, and conditions (a) and (b) in force, 
assume in addition that the family T has the following properties: 

(c’) each F £ 3^ is supported in some closed ball Bp C where the set of radii of 

the balls Bp is bounded away from zero and infinity, 

(d) F G C^{Bp) for any F G T, and sup llVFH^q^ ^ < oo, 

(e) sup ||F||ioo(B^) < oo. 

Then there exist an irrational vector n € Sq, and an unbounded, strictly increasing 
sequence of positive numbers {A^}^^ such that for any Xq G there exists a real-valued 
function vq G C°° (T'^) for which the estimate 


(4.31) 


F(x)uo(AfcM(x, 0) + \kXo)dx - co(uo) 


F{x)dx 


> w(Afc) 


holds for any F G 3" and each integer k'>\, whenever M £ 0{d) and Med = n. 

Proof. For notational convenience we extend all functions F G T to as zero outside 
their supports. Observe that here the finiteness of the supremum of part (b) of Lemma 
4.3 is automatically fulfilled. 

We will start with the case Xq = 0. First carry out the proof of Lemma 4.3 up to 
the definition (4.24). Since now functions from 3“ have compact support, the bound in 
(4.25) can not be obtained directly from integration by parts due to boundary terms 
appearing in the integration. To overcome this technicality we introduce smooth cut-offs. 
Let F G 3" be fixed, and let the closed ball B = B{xo,r) be the support of F. Then by 
(c’) we know that r > cq > 0 for some absolute constant cq. For any A > I/cq we have 
B{xo,r — A“^) C B{xQ,r) and we let ipx : —)• [0,1] be a smooth function such that 
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= 1 on B{xo, r — X), (fx = 0 on \ B{xo, r), and | V</?A(ic)| < ciA for any x G 
where ci is some absolute constant. 

Denoting by Is the characteristic function of the ball B{xo,r), we decompose (4.24) 
into 


I(A;0= [ F(x)Is(x)v9A(x)e2"^"^^^«dx+ / Fix)lBix){l-Mx))e''"'^^"'^^^dx 

:=/i(A;0 + /2(A;0- 

Observe that in /2(A;^) we have integration over B{xo,r) \ B{xo,r — A“^), hence by (e) 
(4.32) |^2(A;0I < l-B(xo,r) \ B{xo,r- A"^)| sup ||i^||Loo(]Rd-i) < CX~^. 

For /i(A; ^) we proceed as in (4.25), however here we will have an additional term coming 
from partial integration, namely the one involving Fd^^px- But observe that all derivatives 
of (fx are supported on B{xo,r) \ B{xo,r — A“^), and hence using the estimate on the 
gradient of (px, along with condition (e) we get 


\F{x)dk(pxix)\dx < ciX\B{xo,r) \ B{xo,r - X 


-1 


)l sup ||F||^o 
Fes' 


i-i) ^ C 


where constants are uniform in F and A. The last bound combined with (4.32) enables 
us to obtain the estimate in (4.25) with possibly different absolute constants. Then, the 
proof of the current lemma for Xq = 0 follows from exactly the same argument in Lemma 
4.3 starting from (4.25) up to the end. 

We now consider the case of any Xq G M'^. Let n be the normal, {Afc} be the sequence, 
and uq be the function for which (4.31) holds with Xq = 0. By constructions of Lemma 
4.3 (Step 3 in particular) and for the case of Xq = 0 there is a strictly increasing sequence 
of integers {im}m=i a set A = C \ {0} satisfying for 

all m G N, such that 


M0) = Y1 

m=l 



'^2TTi^(^)-e g-27rig(’")-6l' 


e G T'^. 


We now slightly adjust the coefficients of vq to handle the effect of the shift. Namely, 
consider the function 

OO 

m=l 





,,-27riA™|(’")-Xog27rj$(™)-0 g27ra,„.j(™)-Xog-27ri^(™)-6» 


By definition, the Fourier coefficient c^{vq) is the complex conjugate of c_g(uo) for any 
^ G Z'^, hence vq is real-valued. It is also clear that vq G C'°°(T'^) and co(uo) = 0. Following 
(4.19) and plugging vq into (4.31) for each integer fc > 1 we get 

J F{x)vo{XkM{x,0)+XkXo)dx = j F{x) cos{2'KXkx'^i'^^^)dx+X:i{Xk)+F. 2 {Xk), 

Rd-l ]Rd-l 


where Si and S2 are defined in analogy with (4.19). Observe that the integral on the 
right-hand side of the above equality is the same as for Xq = 0, and the sums in Si and 
S2 can be estimated exactly as in the case Xq = 0. Indeed, the only difference is that 
coefficients in the sums are multiplied by complex numbers having length 1 (namely the 
exponents involving Xq). This fact will have no effect when taking absolute values of the 
terms in the sums, which is precisely what we do to bound Si and S2. Since the analysis 
is reduced to the case Xq = 0, the proof of the lemma is now complete. □ 
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Proof of Theorem 1.1. We will assume that 
(4.33) ^ oo, as t —)■ oo. 


This assumption results in no loss of generality, for a similar argument as in (4.23). The 
reason for (4.33) is to have slower speed of decay than the error term involved in (4.8) for 
a parameter k = l/(4d). 

Let 1 < 7 < d be fixed from (1.5). From (4.1) and (1.5) we get x*’"' = 0- The latter 
combined with (4.2) implies that for the corresponding boundary layer corrector we have 
= 0 for any n G For 1 < a < d and n G Vn^ solves a uniformly elliptic 

PDF in Q.n with periodic and smooth coefficients, and with boundary data x*’“, hence 
standard elliptic regularity implies that there is a constant Cq independent of n, such that 
|Vynn“(y)| < Cq for any y G dVLn and each 1 < a < d. From this, and the fact that 
Vn^ = 0 it follows that there exists an open subset of the sphere S..,, C such that for 

any n € and any M G 0(d) with Me^ = n, one has 


(4.34) 


1 + ria [a^^x*’"(A(z', 0)) + 0))] 


> 


2 ’ 


for all z' G and any A > 0. Indeed, we simply choose so that each n G S-,, has its 
7 -th component sufficiently close to 1. 

For n G we fix some M„ G 0(d) satisfying M^ed = n and let G°’"'(-,-) be 

the Green’s kernel for the pair M!^). Consider the family of functions 3" := 

{TnlneS'*-!) where we have Fn{x) := d 2 ^dG^’"'{&di (^^jO)) for x G By Lemma 4.2 the 

family T satisfies conditions (a), (b) and (c) of Lemma 4.3. 

Recall that solutions to boundary layer problems are constructed via the reduced 
boundary layer systems of form (2.2), hence we have 

(4.35) vT{y) = vTiMnz) = v^iz) = V^^(N„z', Zd) 


where solves the corresponding problem (2.2), and as is usual = [Nn\n]. In 
particular we have that is Z'^-periodic for any t > 0, and is smooth with respect 

to all its variables. Here one should take into account the snbtlety, that , and hence 
also Vn“, implicitly depend on the matrix M„, but as the choice of M„ is now fixed, we 
may ignore this dependence. 

For n G consider the function 


^n{y) :=l + n„(n-V,x“(y) + aiC(y>0)), y G M'", 

where is fixed from (4.35). Clearly, G C'°°(T'’*). From (4.6) we see that = n-Vy 
which gives the relation between normal derivatives. Now, if y G we get 

(4.36) dtV:{Nz',^) = dz,^T{z,Q) = n-VyV*^^{y). 

From here, the definition of 4;'^ and (4.34), let us show that for any irrational n € E>^ one 
has 


(4.37) |^n(y)|>l/2 for all y G M'^. 

The small nuance, that (4.37) needs the normal to be irrational as compared with (4.34) 
lies in the fact that (4.34) holds on the boundary of 12^, while here we need the entire 
space To see (4.37), observe that for y = Nz' with z' G the lower bound we need 
is due to (4.34) and (4.36). Now, if the normal n is irrational, then {Nz' : z' G is 

everywhere dense in T'’*, which is the unit cell of periodicity of '!'„) hence the continuity 
of completes the proof of (4.37). 

We now apply Lemma 4.3 for the family 3" and modulus of continuity 2oj, and let v € 
be the unit irrational vector and be the increasing sequence given by Lemma 4.3. 
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Thus for v we have (4.37). Next, for a function Fy{x) let vq G C°°(T'^) be the function 
given by Lemma 4.3 for which 


(4.38) 


F^{x)vo{XkM^{x, 0))dx - co(uo) 


Fy{x)dx 


> 2a;(Afc), 


where co(uo) is the 0-th Fourier coefficient and k = 1,2,... . Ellipticity of A implies that 
A{y)v > co|z^P = Co for any y G M'^, with absolute constant cq > 0, hence taking into 
account (4.37) we define 

(4.39) := y e 

and get vq G C°°(T'^). 

Finally, we claim that v, and uq defined by (4.39) satisfy the Theorem. Indeed 

by (4.8) the solution to boundary layer problem with these parameters has the form 

v{y) = [ d 24 G°’''ied,iz', 0 ))vo{XM^{z', 0 ))dz'+ 0{X~^), 

Js.d-1 

where the parameter k in (4.8) is set to l/(4d). The last expression combined with (4.38) 
and (4.33) completes the proof of the Theorem. □ 


5. Application to boundary value homogenization 

This section is devoted to the proof of Theorem 1.3, but before delving into details, 
we sketch the main idea behind the proof. By [2]-[4] we know that the boundary value 
homogenization of type considered in (1.7) is determined by geometric properties of the 
boundary of the reference domain, such as non-vanishing Gaussian curvature [2]-[3], or 
flat pieces with Diophantine normals [4]. Under these conditions one is able to deduce 
effective upper bounds on convergence rates for the homogenization, where the rate will 
be uniform with respect to the boundary data. With these in mind, a suitable candidate 
of domain D for Theorem 1.3 is a C°° domain such that part of its boundary has non¬ 
vanishing Gaussian curvature, while the rest is a piece of a hyperplane. Then relying on 
integral representation of solutions to (1.7) via Poisson kernel, one splits the integral into 
two parts, namely over curved and flat boundaries. The next step is to show that the 
contribution of the curved part has a prescribed rate of decay determined only by the 
embedding of dD into M'^, and hence is invariant under rotations of the domain. This 
step can be fulfilled by adapting the methods of [2]-[3]. For the integral over the flat 
part one shows, using Lemma 4.5, that after a suitable rotation of the domain and an 
appropriate choice of the boundary data it can be made comparatively large. In this 
section we rigorously implement this idea. 

5.1. Preliminary results. We present some technical results which will be used for the 
proof of Theorem 1.3 below. 

Assume we have a bounded domain D C {d > 2) with smooth boundary, and a 
divergence form operator L := —V • A(x)V where coefficient matrix A is defined in D, 
and is strictly elliptic, and smooth. Note, that we do not impose any structural condition 
nor any periodicity assumption on A. Next, we let P(x,y) : D x dD —)> M be the Poisson 
kernel for the operator L in the domain D. Then by Lemma 5.6 we have 

(5.1) \P{x,y)\ < Cp xeD,yedD, 

\x — y\°‘ 

where the constant Cp = Cp{A, D, d). With this notation we have 
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Lemma 5.1. Let 11 he an open and connected subset of dD and Dq D be fixed. Then 

inf / P{x,y)da{y) > 0. 

*S-Do Ju 

Proof. As a trivial observation, before we start the proof, notice that if 11 is the entire 
boundary of D, then the integral in question is identically 1. The general case, however, 
requires some care. The proof is motivated by [3, Lemma 3.1]. We will use integral 
representation of solutions, to get a more precise version of the maximum principle. Fix 
a sequence of smooth functions {gn}’^=i, where gn ■ ^ [0,1] such that gn = I on 11, 

and for any domain 11 C which compactly contains 11, the sequence {gn} uniformly 
converges to 0 outside IT as n —>• oo. We let Un be the solution to Dirichlet problem for L 
in domain D having boundary data gnlgj^- Fix some ^ G IT. Since FI is open in dD there 
exists <5o > 0 such that 


(5.2) 

We get 

(5.3) \unix) - gniOl = 


{yedD-.\y-^\< <5o} C U. 


IdD 


P{x,y)[gniy) - gniOWiv) < 


IdD 


\P{x,y)[gniy) - gn{0]\dcr{y) = 




\P{x,y)[gn{y) - gn{0]\do-iy) < 2Cpd{x)\\g. 


n\\L°° 


'h-€l>'5o 


dcrjy) 
\x - yf 


where we have used (5.2) and the fact that ^(,1 = 1 on IT to pass from the second row of 
(5.3) to the first expression of the last row. We now choose x G D such that |x —< do/2. 
The triangle inequality implies \x — y\ > 5o/2 for all y G dD satisfying \y — f,\ > 6o. Hence, 
the last integral in (5.3) can be estimated as follows 


(5.4) 


da{y) 


\d - 


< 2 ° 


d(y{y) < 


.1 ^d-2 




'\y-^\>So I® y\^ 4|2;-5|><5o \y ?r JSo 

with some positive Cq = Co{d) uniform in x and do, and we have invoked integration 
in spherical coordinates to estimate the surface integral. Without loss of generality we 
assume that constants Co,Cp > 1. We now fix xq G H such that 

do ^ u \ ^ do 


(5.5) 


ko -^1 < Y 


and 


lOCoCp 


< d{xo) < 


ACoCp ■ 


This is always possible, provided Jq > 0 is small enough. Denote D := Dq U {xq}. From 
(5.4) and (5.3) we obtain |m„(xo) — 5n(0l ^ 1/2) hence the triangle inequality implies 

(5.6) |wn(a:o)| = \gnif) + Un{xo) - 5n(0l > \9n{0\ “ kn(a:o) - 5n(0l > 1 - 1/2 = 1/2. 

From (5.6) one has sup^ Un >1/2 for any n G N. From the maximum principle we infer 
that all Un are everywhere non-negative in D, thus applying Moser’s version of Harnack’s 
inequality (see e.g. [14, Theorem 8.20]) we get 

(5.7) 1/2 < sup Un < Cl inf Un < Cl ini Un, n = l,2,..., 

D D D 


where the constant Ci = Ci{d, A, D, D). We have the representation 


Unix) 



Pix,y)gn{y)daiy), 


X G D. 
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Using the construction of gn we pass to the limit in the last integral, getting by (5.7) that 
I P{x,y)da{y) = lim u„(x) > Vx G Dq. 

Jyi n^oo 

The proof is complete. □ 


We will also need a version of the last lemma for a family of Poisson kernels correspond¬ 
ing to a rotated images of a given domain. To fix the ideas, recall that the coefficient matrix 
of (1.7) is defined on some fixed domain X. Let D <£ X he a bounded domain with C°° 
boundary. For a matrix M £ 0{d) define an orthogonal transformation M : —)• by 
Mx = Mx, X G and consider the rotated domain Dm '■= X[D. Obviously Dm is a 
bounded domain, it is also clear that SDm = M(clllo) since M is a diffeomorphism. Next, 
the smoothness of Dm follows directly form the definition of smooth boundary (see e.g. 
Section 6.2 of [14]). We shall take the diameter of D sufficiently small so that Dm <£ X 
for any M G 0{d). Now we get the following version of the previous lemma. 


Corollary 5.2. Let D <Z X be as above and let 11 C dD he open and connected. Fix 
Bq D and for some small constant cq > 0 denote 

O := {M G 0{d) : Bq C Dm and dist{Bo, ODm) > cq}. 


Then 


jnf 

Meo, x&Bq 



PM{x,y)da{y) > 0, 


where Hm = M(n) and Pm is the Poisson kernel for the pair {A,Dm)- 


Proof. For each fixed M G O the infimum is positive in view of Lemma 5.1, and we 
will simply follow the dependence of constants in the proof of Lemma 5.1 on the rotation 
introduced by M. First of all by Lemma 5.6 we have that the constant in (5.1) is indepen¬ 
dent of M therefore (5.5) is uniform with respect to M. Concerning the use of Harnack 
inequality in (5.7) referring to [14, Theorem 8.20] we know that constants for a ball of 
radius R > h depend on dimension of the space, ellipticity bounds of the operator, and 
the radius R. Here we have the same operator for all domains Dm- Finally, the uniform 
distance of Bq from the boundary of Dm for each M G O, along with standard covering 
argument extending the Harnack inequality from balls to arbitrary sets, shows that the 
constant Ci of (5.7) can be chosen independently of M. As the choice of all constants 
in the proof of Lemma 5.1 can be made uniform with respect to M G O, the proof is 
complete. □ 


The next lemma is used in the localization argument of Proposition 5.14 below. 

Lemma 5.3. For ro > 0 let if £ C^{Bro{0)) and assume that V'(O) = |VV'(0)| = 0 and 

HessV'(O) := (9^;3^)l/3=i(0) = diag(ai, ...,ad) £ Md(M), 

where 0 < ai < ... < a^. Then there exist positive constants cq = co(d, [[V'llcs) < = 

ci(rf, IIV'llcs) < ro/ai, and Ki = Ki{d) < 1 < K 2 = K 2 {d, [[V'llca) such that 

(i) Kiai\x -y\ < |VV'(x) - VV'(2/)| < K 2 \x - y\, Vx,y G H(0,ciai), 

(ii) if day is the Kronecker symbol, then for any 1 < a, fi < d one has 

\dly'4^{x) - aaday\ < Vx G H(0,ciai), 

(hi) H(0, coaf) C (VV’)(H(0, cioi)). 
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Proof. We start with part (i). For any x,y £ -Bro(O) by Mean-Value Theorem we have 

d d 

(5.8) |vV'(x) - vv^(2/)| < |9„V’(a;) - 9„V’(y)l < Y - vl 

ol=1 q:,/3=1 

which demonstrates the upper bound of (i). To obtain a lower bound, for 1 < a < d 

d 

and X £ B^o(O) set gaix) = \dl^'ijj{x)\ - Y1 obviously gdO) = > 0. By 

/3=l,/3^a 

C^-smoothness of 'll: we have that each ga has linear modulus of continuity, hence there 
exists a constant ci = ci(d, HV'Hca) such that for all |x| < ci 5 q,( 0 ) we get gaix) > g'a(0)/2. 
From here for any x,y £ (0) Mean-Value Theorem yields 

d d 

(5.9) IVlpix) - V'llj{y)\ > CdY - 9 a'lpiy)\ = Cd ^ |V(5aV’)('ra) • (x - y)\, 

OL=l a=l 

where Tq lies on the segment [x, y] and Cd is a constant depending on dimension. We next 

fix 1 < a < d from \xa — ya\ = max \x^ ~ d/?!) and invoking (5.9) we get 

l</3<(i 

(5.10) |VV’(x) - VV’(y)| > Cd|V((9aV')(rQ,) • (x - y)\ > Cdgaira)\Xa - Val > CrfUilx - y\. 
Combining (5.8) and (5.10) for any x,y £ i?ciai(0) we obtain 

(5.11) Kiai\x -y\< |VV’(x) - VV’(y)| < K2\x - y\, 

with constants K2 = dC 2 (d, ||V'llc 2 ) > 1 and Ki = Ki{d) < 1. This completes the proof of 
part (i). 

Since for each 1 < a, /3 < d has linear modulus of continuity, the claim of (ii) 

follows easily by Mean-Value Theorem with ci sufficiently small. We now proceed to (hi). 

By (5.11) the mapping is invertible in a neighbourhood of the origin, and (iii) is 
simply an effective version of Inverse Mapping Theorem. The desired bound follows from 
the estimate ||(HessV'(0))“^|| x a]”^, C^-smoothness of 'ijj, and [21, Theorem 1.1]^. 

The proof of the lemma is complete. □ 


5.2. A prototypic domain. We introduce a class of domains, call them prototypes, which 
will be used in the proof of Theorem 1.3. Let T be a convex polytope, i.e. a convex 
bounded domain, which is an intersection of a finite number of halfspaces. We assume 
that T C {x £ : Xd < 0} and that 0 £ is an inner point of 87 n {xd = 0}. We fix 

some Ho d 87 n {xd = 0}, a (d — l)-dimensional closed ball centred at 0. Now let Dq C 7 
be a bounded domain having the following properties: 

(PI) Dq is convex with C°° boundary, 

(P2) 8 D 0 n {x £ : Xd = 0} = Ho, 

^Observe that the function V’(x) = ctix? + + o-dX^ manifests that order of the radius of the ball in (iii) 

is generally the best possible. Also, since ip £ , here one may have a direct treatment for (iii) by a 

well-known approach to Inverse Function Theorem. Indeed, set F{x) = V'i/'(x), |a;| < ciui and let j/ € R"* 
be fixed. We need to determine the range of y where the equation F{x) = y has a solution in x from 
B(0,ciai). For this, one may utilize the celebrated method of Newton for finding roots of equations by 
studying the mapping G{x) = x — (VF(0))~^(F(x) — y), where the Jacobian of F is the Hessian of ip. 
Clearly F{x) = j/ iff G{x) = x, i.e. it is enough to figure out when G is a contraction. The latter can be 
resolved easily relying on the G^-smoothness of ip, and determining the range of y when G maps the closed 
ball 5(0, coffli) into itself and has differential of norm less than 1. The details are easy to recover and we 
omit them. 
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(P3) if X G dDo with Xd 7 ^ 0, the Gaussian curvature of dDo at x is strictly positive, 
(P4) we fix some ball Bq lying compactly inside Dq. 

Typically we will embed the whole construction inside a given large domain X. Exis¬ 
tence of Dq satisfying (P1)-(P4) follows directly, as a special case, from an elegant con¬ 
struction due to M. Ghomi in connection with smoothing of convex polytopes, see [13, 
Theorem 1.1], The following picture gives a schematic view of the construction. 



Figure 1. A prototyplc domain Dq obtained as smooth approximation of a polytope. 
Here X is some fixed domain containing 0 € in its interior. Then T is any convex 
polygon sitting inside A n {a; £ < 0}, with non-empty interior, and with part of 

its flat boundary lying on the hyperplane {xd = 0}. We next take a closed flat ball Ho, 
the dashed part on the boundary of T, and invoke [13, Theorem 1.1]. Finally, a ball Bq 
is fixed compactly inside Dq. 


The following notation will be in force throughout the section. Set Tq := dDo and for 
(5 > 0 denote T^ = {x G Pq : dist(x,no) > 5}, where Ho is the {d — l)-dimensional ball 
fixed from (P2) above. Define 

(5.12) k(6) = min min Kq,(x), 

xePji<o<d-i 

where Kq(x) is the a-th principal curvature of Pq at x, and the minimum over P^ exists 
in view of the smoothness of Pq and compactness of P^. In the sequel we assume <5 > 0 is 
small enough so that P^ 7 ^ 0. Due to property (P3) we have 

(5.13) k{6) > 0 and k{6) \ 0 as 5 —)• 0 -|- . 

The next proposition is one of the key ingredients of the proof of Theorem 1.3. 


Proposition 5.4. There exists a modulus of continuity cjq determined by the deeay rate 
of the funetion k{ 6) defined in (5.12) such that for any smooth function P : Pq —>■ any 
g G C'°°(T'^) satisfying fjd 5 = 0, any yo G and any M G 0{d) one has 


(5.14) 


To\no 


P{y) 9 {>^My + yo)da{y) 


< C'a;o(A)||P||ci(ro)ll5'llc‘^(T‘i)> 


VA > 1, 


with a positive eonstant C depending only on dimension d and embedding o/Pq in 


Remark 5.5. One may elaim a decay of the integral in (5.14) relying on [18] for example, 
however without any explicit bounds we have in the current formulation and which we need 
for applications. The proof of this proposition is based on adaptation of methods from 
[2]-[3] both of which work with strictly convex domains, showing that integrals similar to 
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( 5 . 14 ) and involving singular kernel (namely, Poisson’s kernel) decay with some prescribed 
algebraic rate as A —)• 00 . The difference of the current case from [2] -[3] is that on one 
hand here we do not have a singularity introduced by an integration kernel, which gives 
an extra freedom to the entire analysis. On the other hand the strict convexity of the 
hypersurface deteriorates, and we have integration over a hypersurface with boundary; 
both of these factors introduce some technical difficulties which entail somewhat refined 
analysis at certain points. 


Proof of Proposition 5.4. The proof is partitioned into few steps. 

Step 1. Localization. We localize the integral of (5.14) in a neighbourhood of each 
point 2 ; G To of positive curvature. Fix <5 > 0 small. The hypersurface Tq is locally a 
graph of a smooth function, thus there exists r > 0 small such that for any z ^T 25 there 
is an orthogonal transformation iR : —)• satisfying 

(5.15) (iR(ro - 2 )) n ^,.(0) = {(x', V’(a^O) : W\<r}, 

where x' = (xi, ...,Xd-i), is smooth on {x' G : \x'\ < 2r}, '0(0) = |V0(O)| = 0 and 
for the Hessian of if we have 

(5.16) Hess0(O) =diag(ai,..., aa-i) G Md_i(M), 

(5.17) 0 < k{6 ) <oi < 02 < ... < Od-i < Cq. 


Here is the a-th principal curvature of Tq at and the lower bound of (5.17) is due to 
(5.12), while the universal upper bound, as well as the bound |Hess0(x')| < Cq, for any 
\x'\ < 2r, are both due to the smoothness of Tq. Also, r is independent of z and 5, and 
depends on the embedding of hypersurface Tq in M'^. 

By Lemma 5.3 there exist constants Ki < K 2 and cq < ci controlled by dimension d 
and C^-norm of if, and independent of the principal curvatures of Tq, such that 

(a) ATiOilx'l < |V0(x')| < K 2 \x) for all \x'\ < cioi, 

(b) for any l<a,/3<d — 1 setting to be the Kronecker symbol, we get 

\d‘ipif{x') - aa5o,g\ < |x'| < cioi, 

(c) \y'\ < coa^ implies that there exists a unique \x'\ < ciai so that Vif{x') = y'. 


Denote 

(5.18) 


Ls 


Cq 

AdK2 


a 


2 

1 


and consider a family of balls T> = {B{z,^Lfij : z ^ T 2 s}. Clearly covers r 25 , hence by 
Vitali covering lemma there exists a finite collection of disjoint balls Tq = {B{zj, \Lfij : 
j = 1,..., M^} C 23 such that 


Ms 

(5.19) r25 C 1^ B{zj,Ls) =: r25. 

i=i 

As Tq is a graph in L^-neighbourhood of any z G T 2 S get volci-i{B{zj, Ls/5) n r 25 ) 
From this, (5.18) and (5.17) it easily follows that 

(5.20) Ms < < C(k(<5))-2(^-i), 

with an absolute constant C. For 1 < j < Ms set Bj := B{zj,Ls), we now define a smooth 
partition of unity subordinate to these balls. Fix a smooth function T : —)■ [0,1], such 

that 4'(x) = 0, for |x| > 1 and 4'(x) = 1 for |x| < 1. Define 4'5j(x) = 4'(L^^(x — Zj)), 
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Ms \ ^ 

X] '^5j- By construction we have 

i=i ’ J 

Ms 

that ifsj is supported in Bj, ^5,j = 1 on r 25 , and 

i=i 

(5.21) \V(ps,j{x)\ < CL'^^, X G 

where the constant C is independent of j and 6. Now for 1 < j < Ms define 

(5.22) Ij ■= [ P{y) 9 {^My + yo)<fs,j{y)da{y). 

JTo 

Step 2. Reduction to oscillatory integrals. In (5.22) make a change of variables by 
setting y = + Zj, where the orthogonal transformation fR is fixed as in (5.15). We 

next observe that the integral in (5.22) is over B{zj,Ls) n To, and Fq is a graph in the 
L^-neighbourhood of zj. With these in mind we make the change of variable in Ij and 
then pass to volume integration obtaining by so 

(5.23) 

I, = j {P^s,j)izj+^-\z',i;iz')))g {XMzj + XMJl-\z',i;{z')) + y^) il+\Vij{z')\^)P^dz'. 

\z'\<Ls 

For t G M set exp(t) := Next, for ^ G let c^{g) be the ^-th Fourier coefficient 

of g. According to the assumption of the proposition we have co{g) = 0. Using the 
smoothness of g and expanding it into Fourier series we get 

(5.24) g {XMzj + XMJi-\z^,'ilj{z^)) + yo) = 

^ cg( 5 r)exp (AC • Mzj + C • 2/o) exp [A3?M^C ' (z’, 'ipiz'))] , 
«6ZA{0} 

where we have also used the orthogonality of M and fR. For C G Z*^ set fRM^C ■= g ■= 
\g\{n',nd) with {n',nd) G By orthogonality of M and fR we have |r/| = |C|. Next, 

define 

F{z') =n' ■ z' + ndi^iz')^ 

and 

<h,(z') = {P^s,j){zj + 'R-\z',^{z'mi + \Vi,{z')\^)P\ 

From the definition of the cut-off <ps^j we have 

(5.25) \\^,\\cr<C\\P\\cMj\ 

uniformly for all 1 < j < Ms with an absolute constant C . With these notation, from 
(5.23) and (5.24) we get 

(5.26) Ij= ^ c^{g)ex.^{Xi-Mzj 

«6ZA{0} 

where 

(5.27) = f 4>j(2;')exp [A|C|F’(2:')] dz'. 

J\z'\<Ls 

Step 3. Decay of Ij. We analyse the decay of each in two distinct cases. 

Case 1. \n'\ > coa1f2. 


then clearly supp^'^j C Bj. Now let tpsj '■= 
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Fix l<Q;<d—Iso that ItIqI = maxi<^<rf_i |n^|, clearly \na\ > \n'\/d. From this, 
definition of F, (5.18) and assertion (a) of Step 1, on the support of we have 

(5.28) \daF{z')\ = \na + nddaip{z')\ > |n„| - \dai){z')\ > 

^ - K^Ls >^- = Z^a\. 

2d 2d 4dK2 4d ' 

Integrating by parts in /j g in the a-th coordinate, and then employing (5.28) and (5.25) 
we get 


14^1 < c{x\m 


-1 


l\z'\<Ls 


On 


d„F 




d2'<C(A|C|)-4f44||P||ci4- , 


with an absolute constant independent of j. From here and (5.18) we have 
(5.29) 14^1 ^^(AieD-illPlIci^-®. 


Case 2. \n'\ < coaf/2. 

Since |(n',nrf)| = 1 and cq is small we have \nd\ > 1/2 and hence 



< coaf. By (c) 


there exists a unique Xq G P(0,ciOi) such that VV'(xq) = —hence VP(xq) = 0. 
Observe that Xg is not necessarily from the support of ‘hj. For 1 < a < d — 1 consider the 
cone 


clearly U^~\Sq, = Now for fixed 1 < a < d—1 take z’ G Qa such that |xq+z'| < cioi. 

Using estimate (b) of Step 1, the facts that \nd\ > 1/2 and VP(xq) = 0, and invoking 
Mean-Value Theorem, for some r on the segment [xg, Xg + F] we get 


(5.30) 


dF , ; 


dF , , 

dF , ,, 


5P A , , , 


— 



— 



> 


d-l 


\nd\ I \daai^{T)Za\ - Yj | > Cai\z'\, 

0=1, 0¥^a 


with an absolute constant C > 0. Since the cones {Cq} cover for each 1 < a < d — 1 

there exists Wq supported in Cq, smooth away from the origin and homogeneous of degree 
0 such that 

d-l 

Y^a{z') = l, Vz' / 0. 

a.=l 

Observe that since each Wq is homogeneous of degree 0, for all 1 < a < d — 1 and non-zero 
z' G near the origin we have 

(5.31) |9Qa;Q(2:')| < C|^, 

with an absolute constant C. 

Now fix a non-negative function h G C°°(M‘^“^) satisfying h{x') = 0 for \x'\ > 2 and 
h{x') = 1 for \x'\ < 1. Setting x' = z' — Xg form (5.27) we get 


ki = 


4 (xo -b x') exp [A|^|P(xo -b x')] dx' := 


where 

if} = [ /i(A^/^x')<hj(xo -b x') exp [A|^|P(xo -b x')] dx. 

dRd-i 
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From the definition of h we have 
(5.32) 

The second part we decomposed as ^ where 

a=l 

= [ uJa{x'){l - h{X^^^x'))^j(xQ + x') exp [A|^|F(xo + x')] dx'. 

jRd-l 

We now invoke partial integration in a-th coordinate, and with the aid of estimates (5.30), 
(5.31) and (5.25) we get 


r(2),a 


< 

rs_/ 



\x'q+x'\<Ls 


1 

x' 


^(1 - h{\^/\'))Lf\\P\\c^ + \^/\dah){\^/\')\\P\\L^ 


dx'+ 


(A|^|) 


-1 


/ ^(l-^(A^/V))||P||Loodx'. 


\x'q+x'\<Ls 


We now use the definition of h, and that of Ls given by (5.18), and employ integration in 
spherical coordinates by so appearing to 


(5.33) 


(2),a ^ (A|g|) 


-1 


hi 


ai 


:aV2l-| 


P\\c^ + 


(Aiei)-' 


-XP^\\P\ 


with an absolute constant® C. By (5.32) and (5.33), along with the definition (5.18) we 
get 


(5.34) 


Ij,^\ < CX-P‘^\\P\\cia^^. 


Step 4. Final estimates. We now put everything together. By [2, Lemma 2.3] we have 

( 

(5.35) \cd9)\<C Y. l|V“5lli2(T^) 

CeZ'^\{0} yaez^, \a\=d 

where V" = 9“^ o ... o d'^’^, jaj = jaij + ... + ja^j, and constant C = C{d). 

Now let 1(A) be the integral in (5.14). From (5.22) and (5.19) we get 



(5.36) 


Ms 

n\) = Y.ii + 


'(ro\no)\r2i 


Piy)9{>^My + yo)da{y). 


The definition of F 25 infers volrf_i((Fo \ Flo) \ 1 ^ 25 ) < 6, hence by (5.36) it follows that 


Ms 


(5.37) 


i.^(A)i ^ i.^ji+’^ii.piiL°°iigiiL°°. 
1=1 


®One may get more precise decay rate in A depending on dimension, cf. [3, p. 76], however the crude 
estimates we have here are enough for our purpose. 
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We thus get 


Ij\ < ^|cg(5r)||4g| (by (5.26)) 


< 

r\j 




-d9)\ 


(Aiei)-'IIP 


I „‘id-8 I \ 

1 (^ 10 ]^ + A 


-l/2||p| 


^-3 

C-iOi 


(by (5.29) and (5.34)) 


< X~^^‘^{K{6))~^\\P\\ci\\9\\cd (by (5.35) and (5.17)), 


where the constants are uniform in 1 < j < Ms and 5 > 0. Using this bound on Ij along 
with estimate (5.20) on Ms, from (5.37) we get 

(5.38) |/(A)| < [A-i/2(^(5))-2U+i)+<5]||p||^i||5||^,, 

for any <5 > 0 small and any A > 1. It is left to optimize the last inequality in <5, for which 
consider the function f{6) := in the interval (0, (5o) where Jq > 0 is small. 

It follows from definition of k{5) in (5.12) that / is continuous and strictly increasing in 
{0,6o), and converges to 0 as 5 —)■ 0+. Hence for each A > 0 large enough there is a unique 
0 < (5 = (5(A) < (5o such that A“^/^ = /((5(A)). Define u;o(A) := (5(A) and observe that 
wo(A) = /“HA”^/^), where /“^ is the inverse of /. It readily follows from the mentioned 
properties of / that uq is a modulus of continuity. Finally, for given A > 0 large, applying 
inequality (5.38) with 5 = (5(A) we get |/(A)| < a;o(A)||P||(^i||( 7 ||(^d completing the proof of 
the proposition. □ 


Proof of Theorem 1.3. Recall that the coefficient matrix of (1.7) is defined on 
some fixed domain X. Take any xq in the interior of X and consider the translated 
domain X — xq. Since the origin lies in the interior oi X — xq we fix Dq d X — xq, any 
prototypic domain constructed in subsection 5.2. Following the notation of this section, 
we let Fq := dD^ and by Hq we denote the flat portion of Fq, which by (P2) is a (d — 1)- 
dimensional closed ball. Finally, we fix a ball Bq from property (P4) formulated above. 

For a matrix M G 0(d) let bounded domain Dm be the rotated image of Dq by M as 
defined in subsection 5.1. Set Vm = QDm, as we have already discussed Fm is smooth 
and we have Vm = 3V[(Fo). Since M is orthogonal, it follows that Hm '■= MHo C SDm 
is a (d — l)-dimensional ball lying in a hyperplane passing through the origin and having 
normal equal to Med- Finally, we fix a small closed neighbourhood Sq C of es with 
non-empty interior, such that for any n G 5o if M G 0(d) with Med = n, then the 
corresponding rotated domain Dm lies in X — xq and compactly contains the ball Bq. 

Fix a function g G C'“(T'^) with the property fjd 9 = 0, and a unit vector n G S'o along 
with a matrix M G 0(d) satisfying Med = n. For e > 0 let be the smooth solution to 

(5.39) — V • A{x)VUe{x) = 0 in Dm + xq and Us{x) = g{x/e) on Vm + a^o- 

Since g has mean zero, by (1.8) we get that uq, the homogenized solution corresponding 
to (5.39), is identically zero. Now let Pm ’■ {Dm + xq) x (Fm + xq) —)• M be the Poisson 
kernel for the pair {A, Dm + xq). For x G Dm + xq we have 

(5.40) Ue{x) = / PM{x,y)g{y/e)daM{y) = / + := h + h- 

9Um+X 0 9rM\^M+^0 

Observe that F^ \ H^ is the part of the boundary of Dm with non-vanishing (positive) 
Gaussian curvature, while VIm is a flat ball (notice that principal curvatures are orthogonal 
invariants, hence do not change after rotation, see for example [14, Section 14.6]). The 
aim is to show that I 2 has some prescribed decay rate in e independently of M, while the 
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decay of Ii can be made as slow as one wish. Setting A = 1/e and then translating the 
origin to xq and rotating the coordinate system by we get 

h = h{^]x)= / Pm{x,Mz + xo)g{\Mz + \xo)daQ{z). 

Jro\no 

For each fixed x £ Bq + xq hy (5.50) the function Pm{x,- + xq) is (7°° on Vm, hence 
applying Proposition 5.4 we obtain 

\I 2 \ < Coujo{X)\\Pm{x, ■ + a;o)||ci(rM)ll5'llc‘^(T‘^)’ x £ Bq + xq, 

with an absolute constant® Co- By Lemma 5.6 we have \\Pm{x,- + xo)||ci(rM) — 
uniformly for x £ Bq + xo and matrix M as above. The latter implies that 

(5.41) |l 2 (A;x)| < C'oa;o(A)||c/||c.<i(Td), x £ Bq + xq. 

For the flat part of the integral we denote by Im the characteristic function of 11^ and 
extend Pm{x, •) as zero outside 11^ + xq- With these notation we get 


(5.42) 


/i=/i(A;x)= / PM{x,y + xo)g{Xy + Xxo)daM{y) = 
JTIm 

/ PAiix, y + xo)lMiy)giXy + Xxo)da{y) = 


' y-n={) 


'Zd=0 


Pm{x, Mz + xo)lMiMz)g{XMz + Xxo)dao{z) = 


/ Pm (x, M (/, 0) + Xo) Im {M (z', 0)) g (AM(z', 0) + Axq) dz'. 

jRd-l 

Consider the following 2-parameter family of functions 


(5.43) S':= {Fx,Miz')'■= PMix,M{z ,0) + xo)lMiM{z ,0)), z'G where 

X £ Bo + Xo, M £ 0(d) with Mcd = n and n £ S'o}. 

Let us see that the family 3" satisfies all conditions of Lemma 4.5. Indeed, the assumption 
(c’) of the lemma is simply due to the construction of domain Do and orthogonality of 
matrices M. Next, (c’) combined with the estimate (5.51) implies assumption (b) of 
Lemma 4.3. Items (d) and (e) of Lemma 4.5 are due to (5.52) and (5.51) correspondingly. 
It is left to verify the assumption with the lower bound on the integrals in part (a) of 
Lemma 4.3. The latter is due to Corollary 5.2. 

We now apply Lemma 4.5 to the family T for modulus of continuity A i—>■ a;(A) -|- 
ujq (A), where w is the one fixed in the formulation of the theorem, and cuq is determined 
from (5.41). By doing so we get an irrational normal n £ So, a function g £ C'°°(T'^) with 
the property g = 0, and a strictly increasing sequence of positive numbers 
such that for any M £ 0(d) satisfying Med = n from (5.42) we have 

(5.44) \Ii(Xk-,x)\>u}(Xk)+ujl^'^(Xk), x£Bo + xo, A: = 1 , 2 ,.... 

Since wq decays at infinity, combining (5.44) and (5.41) and setting = 1/Xk, from 
the representation (5.40) for any x G Bq -|- xq we get 

(5.45) 

1 /2 

\u£^:(x)\ > \Ii(£k-,x)\ - |l 2 (efc;x)| > uj(Xk) + ujq i^k) - Cou}o(Xk)\\g\\cd{T‘i) > ^(l/efc), 

if /c > 1 large enough. Finally, as a domain D of the Theorem we take Dm + xq where 
Dm is any domain having n as the normal of its flat boundary, where n is obtained from 

^Observe that the constant in Proposition 5.4 is independent of the shift and hence having Axq in g still 
results in a constant independent of A. 



32 


HAYK ALEKSANYAN 


Lemma 4.5 (note that D is defined modulo 0{d — 1)), and as D' we take Bq + xq, where 
Bq is the ball fixed above. Since fjd 9 = 0, the homogenized problem (1.8) has only a 
trivial solution. Thus (5.45) completes the proof of part a) of the theorem. 

We now prove part b). Observe that an argument similar to what we had for (5.41) 
shows that x) —>■ 0 as A —>■ oo for any x G Z1 = Dm + a:o- Hence, taking into account 
decomposition (5.40), it is enough to show that Ii(A;x) —)> 0 as A —>■ oo for all x G D. 
We may assume without loss of generality, by passing to a subset of Sq if necessary, that 
Ud 7 ^ 0, where the unit vector n is fixed from part a). Now take any subset H of Bm + xq 
such that the projection of H onto x {0} is a (d — l)-dimensional rectangle, which 

we will denote by [ai,5i] x ... x [arf_i,5rf_i]. To show that /i(A;x) decays it is enough to 
prove that the integral 


(5.46) 


/(A;x) 



PM{x,y)g{\y)daM{y) 


converges to 0 as A —)• oo for any x € D. The latter follows easily. Indeed, since / 0 
for y G n we have yd = —^(niyi + ... + Ud-iyd-i) hence passing to volume integration in 

(5.46) and expanding g into Fourier series, for each x € D we get 

(5.47) 






d-l 

PM{x,yi,...,yd) 

k=l 


exp 



XVk 

'^d ) 


dyi...dyd-i, 


where yd = — ;^(ui?/i + ... + Ud-iyd-i), and as before exp(t) = t G M. Since n is 

irrational, it is easy to see that at least one of the exponentials in the last expression is non¬ 
trivial, i.e. for any non-zero G there isl<A;<d—1 such that — xwCd 7^ 0. Using 
this and the smoothness of Pm{x, ■) which is due to Lemma 5.6, in (5.47) for each non-zero 
^ we integrate by parts with respect to the corresponding k-th. coordinate, obtaining by 
so that each integral in (5.47) decays as A —)• oo. Observe as well that for each fixed 
X, integrals in (5.47) are uniformly bounded with respect to A and ^ in view of (5.51). 
Finally, the series of Fourier coefficients of g converges absolutely due to the smoothness of 
g. This, coupled with uniform boundedness and decay of integrals in (5.47), easily implies 
that /(A; x) —>• 0 as A — )■ oo for each fixed x G H, completing the proof of part b) of the 
theorem. 

The theorem is proved. □ 


5.3. Concluding comments. The reader may have observed that the approach we have 
here has a potential to work for homogenization problems when solutions to the underlying 
PDF admit integral representation with some nice control for representation kernels. For 
example, one should be able to treat the periodic homogenization of Neumann boundary 
data with the methods developed in this note. It also seems plausible, possibly with some 
more work, that one can study homogenization of almost-periodic boundary data as well 
using a similar analysis. 

The reason we can not readily allow (periodic) oscillations in the operator of the prob¬ 
lem (1.7) is due to the absence of the necessary control over the Poisson kernel correspond¬ 
ing to oscillating operator. In particular, we do not have uniform (with respect to e > 0) 
bounds on C^-norms similar to those in Lemma 5.6 at our disposal. One possible detour 
of this obstacle is the use of results on homogenization of Poisson kernel for the e-problem 
considered in [16], to reduce matters to a fixed operator, where the analysis of the current 
paper can be utilized (cf. the proof of Theorem 1.7 of [3]). 
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Appendix 

We give some basic estimates on Green’s and Poisson’s kernels associated with diver¬ 
gence type elliptic operators in bounded domains. The estimates we will need are standard 
and well-known, but in view of the absence of an explicit reference we include the proofs 
here. 

Following subsection 5.1 for a matrix M G 0(d) and domain D C M '’* by Dm we denote 
the rotated image of D by the matrix M. We saw already that Dm is a bounded domain 
with (7°° boundary, moreover ODm = 'M{dD), where the orthogonal transformation M is 
defined through Mx = Mx for x G M'’*. Next, we let the coefficient matrix A and domain 
X be as in the formulation of Theorem 1.3. Finally, we fix a bounded domain D with C°° 
boundary, and such that for any M G 0(d) the domain Dm lies compactly inside X. 

With these preliminary setup we now consider the operator L := —V • A(x)V, x G A, 
and let GMix,y) be the Green’s kernel for L in domain Dm, i-e. for each fixed y G Dm, 
GM{-,y) solves 

-Xx ■ A{x)XxGMix,y) = 6{x - y), x G Dm and GM{x,y) = Q, x G BDm, 
in a sense of distributions, where 5 is the Dirac symbol. We have the following properties. 

(a) Set G^(x,y) := G'j^{y,x), then is the Green’s kernel for the formal adjoint 
to T in domain Dm, i-e. a divergence type operator with coefficient matrix equal 
to 

(b) For any multi-index a G we have 

(5.48) \X^GMix,y)\ < Ga\x - for d-|- |r?T,| > 2, 

where Go, is independent of M. 

The existence and uniqueness of Green’s kernel, as well as property (a) and estimate of 
(b) for each fixed M are proved in [9], which is one of the several papers discussing some 
basic questions around the Green’s kernels for divergence type elliptic operators. The only 
thing that needs some clarification is the choice of the constant in (b) independently of 
M. To see it, we proceed as follows. 

Assume G is the Green’s kernel for the pair [A,D) and let M G 0{d) be fixed. Make 
a change a variables by setting Mx = x and My = y where x,y G D and x,y G Dm- 
By orthogonality of M we have which combined with the uniqueness of the 

Green’s kernel, easily implies that the function 

G{x,y) ■.= G{M^x,M^^ 

is the Green’s kernel for the pair {MA{M'^-)M'^,Dm)- This shows that if Gm is the 
Green’s kernel for the pair [M'^A{M-)M,D), then 

(5.49) Gm{x, y) := Gm{M^x, M^^ 

produces the Green’s kernel for the pair {A, Dm)- In view of the orthogonality of M 
the ellipticity constants as well as bounds on norms, for any A; > 0, of the coefficient 
matrix M"^A{M-)M can be chosen independently of M, hence (5.49) and the proof of [9] 
illustrate uniformity of the constants in (b) with respect to M. 

Now for X G Dm and y G ODm we let PM{x,y) be the Poisson kernel for the pair 
(A, Dm)- One may easily conclude using the divergence theorem that 

(5.50) PM{x,y) =-rF{y)A^{y)XyGM{x,y), whem x G Dm, y G SDm, 
where n{y) is the unit outward normal to ODm at the point y. 
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Lemma 5.6. Let the domain D be as above and M G 0(d) be any. Keeping the above 
notation and assumptions in force we have 

(5.51) \PM{x,y)\ < C \/x e Dm, yy ^ ODm, 

\x — y\°- 

where d{x) is the distance of x from the boundary of D, and 

(5.52) \VyPM{x,y)\<C\x-y\~‘^, Vx € Dm, yy € ODm- 
with constants independent of M. 


Proof. Let us first observe that the estimates on derivatives of P follow trivially from 
the representation (5.50), estimate (5.48), the symmetry property of Green’s matrix (a) 
formulated above, combined with smoothness of the coefficients A and domain Dm- 

We now proceed to the proof of the estimate with distance, for which we will rely on 
uniformity of constants in (5.48) with respect to M. First consider the case of d > 3. We 
now show that for x,y € D with x ^ y one has 


(5.53) 


\GMix,y)\<C 


d{x) 

X — y\^~^ 


By (5.48) we have \Gm{x, y)\ < G\x — y\^ hence (5.53) is trivial if d{x) > ^\x — y\, thus 
we will assume that d{x) < ^\x — y\. Now fix x G dD such that d{x) = \x — x\. Since 
Gm vanishes on the boundary of Dm with respect to both variables, using Mean-Value 
Theorem and estimate (5.48) for G*j^ we get 


(5.54) \Gm{x, y)\ = \Gm{x, y) - Gm{x, y)\ < {V^GMix, 


X — x\ 


< G- 


d{x) 


x-y\ 


d-i' 


Here, for estimating the derivative of Gm we have used symmetry relation (a) above. As 
d{x) < ||x — ?/|, and x is on the segment [x,x], the triangle inequality implies 

\x — y\'>\x — y\ — \x — x\>\x — y\ — \x — x\ > -|x — y|, 

O 

which combined with (5.54) gives (5.53). 

Now fix xo,?/o and let r := \xo — yo\ > 0. Consider Gm{z) '■= GM{xo,rz + xq) in a 
scaled and shifted domain D^^m '■= '<'~^[Dm — xq). We have that Gm is a solution to the 
adjoint operator in D^-^m H (H(0,3) \ B{0, 1/3)), hence in view of the smoothness of the 
coefficients, from standard elliptic regularity estimates we obtain 


I'^zGMiz)] — \rKyGM{xo,rz + xo)\ < C'||G'M||L°o(D^,Mn(B(o,3)\s(o,i/3)))) 

for all z G Dr^M H (H(0, 2) \ -6(0,1/2)). Note that the regularity estimates we have used 
are uniform in M due to (5.49) and orthogonality of M. The last inequality combined 
with (5.53) implies 

(5.55) \VyGM{xo,yo)\ < 

Now the estimate (5.51) follows from the last inequality and representation (5.50). 

It is left to consider the 2-dimensional case, which can be handled precisely as in [7, 
Lemma 21] thus we will omit the details. 

The lemma is proved. □ 


Acknowledgements. I thank Christophe Prange for an important conversation back in Fall of 2014 which 
has eventually led to the problem considered in Section 3. Part of the work was done during my visit to 
Institut Mittag-Leffler for the term “Homogenization and Random Phenomenon”. The warm hospitality 
and support of The Institute is acknowledged with gratitude. This article is partially based on my PhD 








SLOW CONVERGENCE PHENOMENON 


35 


thesis completed at The University of Edinburgh in 2015, and I wish to thank my thesis supervisor Aram 
Karakhanyan for his encouragement, as well as advice. I am also grateful to the anonymous referees for 
thorough reading of the manuscript and providing valuable suggestions, as well as corrections which have 
certainly helped to improve the quality of the presentation. 

References 

[11 Aleksanyan, H.: Regularity of boundary data in periodic homogenization of elliptic systems in layered 
media. arXiv:1409.7344v3 (2014) 

[2] Aleksanyan, H., Shahgholian, H., Sjolin, P.: Applications of Fourier analysis in homogenization of 
Dirichlet problem I. Pointwise estimates. J. Differential Equations 254(6), 2626-2637 (2013) 

[3] Aleksanyan, H., Shahgholian, H., Sjolin, P.: Applications of Fourier analysis in homogenization of 
Dirichlet problem. estimates. Arch. Ration. Mech. Anal. (ARMA) 215(1), 65-87 (2015) 

[4] Aleksanyan, H., Shahgholian, H., Sjolin, P.: Applications of Fourier analysis in homogenization of 
Dirichlet problem III. Polygonal domains. J. Fourier Anal. Appl. 20(3), 524-546 (2014) 

[5] Allaire, G., Amar, M.: Boundary layer tails in periodic homogenization. ESAIM, Control, Optim. 
Calc. Var. 4, 209-243 (1999) 

[6] Armstrong, S., Shen, Z.: Lipschitz estimates in almost-periodic homogenization. Comm. Pure Appl. 
Math, (to appear) arXiv:1409.2094 

[7] Avellaneda, M., Lin, F.: Compactness methods in the theory of homogenization. Comm. Pure Appl. 
Math. 40(6), 803-847 (1987) 

[8] Bondarenko, A., Bouchitte, G., Mascarenhas, L., Mahadevan, R.: Rate of convergence for correctors 
in almost periodic homogenization. Discrete Contin. Dyn. Syst. 13(2), 503-514 (2005) 

[9] Dolzmann, G., Muller, S.: Estimates for Green’s matrices of elliptic systems by theory. Manuscripta 
Math. 88(1), 261-273 (1995) 

[10] Dong, H.; Kim, S.: Green’s matrices for second order elliptic systems with measurable coefficients in 
two dimensional domains. Trans. Amer. Math. Soc. 361, 3303-3323 (2009) 

[11] Gerard-Varet, D., Masmoudi, N.: Homogenization in polygonal domains. J. Eur. Math. Soc. (JEMS) 
13(5), 1477-1503 (2011) 

[12] Gerard-Varet, D., Masmoudi, N.: Homogenization and boundary layers. Acta Math. 209, 133-178 

( 2012 ) 

[13] Ghomi, M.: Optimal smoothing for convex polytopes. Bull. London Math. Soc. 36(04), 483-492 (2004) 

[14] Gilbarg D., Trudinger N.: Elliptic partial differential equations of second order. Second edition. 
Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences], 
224. Springer-Verlag, Berlin, 1983. xiii-|-513 pp. 

[15] Hofmann, S.; Kim, S.: The Green function estimates for strongly elliptic systems of second order. 
Manuscripta Math. 124, 139-172 (2007) 

[16] Kenig, C., Lin, F., Shen, Z.: Periodic homogenization of Green and Neumann functions. Commun. 
Pure Appl. Math. 67(8), 1219-1262 (2014) 

[17] Kenig, C., Lin, F., Shen, Z.: Convergence rates in for elliptic homogenization problems. Arch. 
Ration. Mech. Anal. (ARMA) 203(3) 1009-1036 (2012) 

[18] Lee, K.-A., Shahgholian, H.: Homogenization of the boundary value for the Dirichlet problem. 
arXiv:1201.6683 (2012) 

[19] Prange, C.: Asymptotic analysis of boundary layer correctors in periodic homogenization. SIAM J. 
Math. Anal., 45(1), 345-387 (2012) 

[20] Schulze, B.-W., Wildenhain G.; Methoden der Potentialtheorie fiir elliptische Differentialgleichnungen 
beliebiger Ordnung; Lehrbiicher und Monographien aus dem Gebiete der Exakten Wissenschaften: 
Mathematische Reihe 60, Birkauser-Verlag, Basel (1977) 

[21] Wang, X.: Convergence of Newton’s method and inverse function theorem in Banach space. Math. 
Comp. 68(225) 169-186 (1999) 

School of Mathematics, The University of Edinburgh, JCMB The King’s Buildings, Peter 
Guthrie Tait Road, Edinburgh EH9 3FD 

Current address: Department of Mathematics, KTH Royal Institute of Technology, SE-100 44 Stock¬ 
holm, Sweden 

E-mail address: hayk.aleksanyan@gmail.com 



